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Drei Beispiele aus der Elektroservotechnik 
Zusammenfassender Bericht 


Von EDUARD GERECKE, Ziirich?!) 


Einleitung 


: In einer fritheren Arbeit wurde in dieser Zeitschrift die theoretische Seite 
der Berechnung von Regulierungen?) behandelt. Die drei folgenden elemen- 
_taren Beispiele sind charakteristisch fiir die Elektroservotechnik. Sie zeigen. 
- zanachst das statische Verhalten des geschlossenen Regelkreises in graphischer 
_ Darstellung. Aus den Kennlinien der einzelnen Teile, wie Generatoren, Motoren, 
_ Apparate und Verstarker, wird die aussere Kennlinie des geschlossenen Regel- 
_ kreises mit Proportionalregler konstruiert. 

Heute interessiert sich die Technik immer mehr fiir das dynamische Ver- 
-halten gegengekoppelter Systeme, und zwar nicht nur hinsichtlich der Stabili- 
tat, sondern vornehmlich wegen der Méglichkeit des Baues von genau und 
rasch arbeitenden, also mit geringen Zeitverzdgerungen behafteten Servo- 

mechanismen. 
Es wird gezeigt, wie in vielen Fallen durch das Schliessen des Regelkreises 
die Zeitkonstanten reduziert werden und welches offensichtlich der Grund 
hiefiir ist. An dem Zahlenbeispiel eines grossen Walzmotors werden die ver- 
 schiedenen Stabilitatskriterien besprochen. Die Arbeit bedient sich vornehm- 
lich graphischer Methoden, da diese dem Ingenieur bei der Projektierung und 
Dimensionierung der Einzelteile wertvolle Dienste leisten. Seine Auigabe ist 
es, mit konstruktivem Geschick technisch einfache und wirtschaftlich tragbare 

| Lésungen zu finden. Die beschriebenen Methoden lassen auch die Behandlung 
komplizierter Probleme zu, wie zum Beispiel das dynamische Verhalten von 

_ grossen Drehstromgeneratoren mit den zugehérigen Turbinen und Fernleitun- 
gen. Treten jedoch nichtlineare Glieder auf, so sind solche Systeme heute einer 
mathematischen Behandlung noch nicht zuganglich. 

Es wurde durchgehend die vom Schweizerischen Elektrotechnischen Verein 
| empfohlene Schreibweise verwendet?). 


1) Institut fiir Allgemeine Elektrotechnik der ETH. 

2) Hernricu H. WEBER, Methodik der Berechnung von Regulierungen, ZAMP 4, 233 (1953). 

3) Regeln und Lettsdtze fir Buchstabensymbole und Zeichen, Publikation Nr. 192df, 2. Auflage 
_ des Schweizerischen Elektrotechnischen Vereins (SEV) (Ziirich 1953). 
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A. Fremderregte Gleichstrommaschine 


Als erstes Beispiel soll die graphische Behandlung der Regelung der Klem- 
menspannung U eines mit konstanter Drehzahl m laufenden fremderregten 
Gleichstromgenerators nach Figur 1 besprochen werden. Figur 2 zeigt die 
Verhiltnisse fiir den offenen und geschlossenen Regelkreis im stationaren Zu- 
stand. Die dusseren Kennlinien AB des offenen Kreises (k = 0) sind bei Ver- 
nachlassigung der Ankerriickwirkung parallele Gerade mit einer Neigung ent- 
sprechend dem «Quellenwiderstand» R, des Ankerkreises. Im geschlossenen 
Regelkreis wird die Klemmenspannung U iiber ein Potentiometer mit dem 
Teilverhaltnis k dem Verstarker mit der Ubertragungsfunktion K, g(p) in Ge- 
genschaltung zur frei wahlbaren Vergleichsspannung U, zugefiihrt. Im Leer- 
lauf (I = 0) fihrt das graphisch in Figur 2 vom Punkt A zum Punkt A, auf der 
«Potentiometergeraden K,» und dann zum Punkt A, auf der statischen « Ver- 
_starkerkennlinie K,», bei welcher vertikal dessen Eingangsspannung E = U,—kU 
und horizontal dessen Ausgangsstrom J, eingetragen sind. Zum Erregerstrom J, 
des Generators kann die zugehérige Quellenspannung U, der «Leerlaufkenn- 
linie K,» im Punkt A, entnommen werden. Der Polygonzug AA,4,4,A stellt 
daher die elektrischen Verhaltnisse im geschlossenen Regelkreis bei Leerlauf 
dar. Bei Belastung im offenen Regelkreis durch den Strom J sinkt die Klem- 
menspannung U langs der Geraden AB laut 


Ue U,— Rol : (1) 


um den Betrag AB) = AU = R,I ab. 

Wie verandert sich nun diese Kennlinie im geschlossenen Regelkreis ? Liasst 
man den Spannungsabfall AU zu, so sinkt die Gegenspannung k U am Verstar- 
ker um (k AU) (Punkte B, und B,). Demzufolge tritt ein um AJ, grésserer 
Erregerstrom auf, der eine gréssere Quellenspannung (U, + 4U,) entsprechend 
den Punkten B; und B, hervorruft. Damit stellt sich ein viel grésserer Bela- 
stungsstrom J* (Punkt B*) ein, wobei unter Vernachlassigung der Ankerriick- 
witkung B,B* || AB ist. 

Durch das Schliessen des Regelkreises geht also die urspriingliche Kenn- 
linie AB in die viel flachere Gerade AB* iiber, wie wenn der Quellenwider- 
stand R, kleiner geworden ware. Die entsprechende Konstruktion fiir den 
hohergelegten Leerlaufpunkt C ergibt zufolge der Kriimmung der Kennlinien 
kK, und K; eine im Vergleich zu A B* steilere Kennlinie CD*. Auf ein Absinken 
der Klemmenspannung U um AB, = AU reagiert also der Regelkreis durch eine 
Erhohung der Quellenspannung U, um AU,; Figur 3 ist typisch fiir den «Pro- 
portionalregler». Die Ubertraeungsfunktion. K, g(f) des Verstarkers habe fiir 
p = 0 den Wert K,. Der Erregerstrom J, in der Feldwicklung mit dem Wider- 
stand R,, der Induktivitat L, und der Zeitkonstanten 7, = L,/R, wachst infolge 


a 
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Figur 1 
Gleichstromgenerator mit gegengekoppeltem Regelkreis. 


Figur 2 
Aussere Kennlinien des Generators bei offenem und geschlossenem Regelkreis. 


Figur 3 
Fehler und Quellenspannungsanderung beim Proportionalregler. 
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der Spannungsabsenkung AU um 


Pee 
[=k AU = = AU , : 


was unter Beriicksichtigung der Ubertragungsfunktion C des Generators eine | 
Erhéhung der Quellenspannung um : 


AU, = AB,= C AI, = —— 4U=K AU 
hervorruft. Nach Figur 2 ist damit: 


BB, =U + AUS 4K) AU, BAA eee ee 


VERE: 
fae ae eee 


1+ Kk’ Re 2) 


Der geschlossene Regelkreis reduziert also den Quellenwiderstand R, auf einen 
(1 + K)-mal kleineren Wert R>. 

Die Kennlinien des geregelten Kreises setzen sich nach Figur 2 aus einem 
links der Grenzgeraden G,G, gelegenen geraden Teil und aus einem durch die 
Kriimmung der Kennlinien K, und K, bedingten gekriimmten Teil zusammen. 
Sie sind zueinander kongruent und lassen sich durch Translation langs G,G, 
zur Deckung bringen. Statt der Konstanten K und C ist im gekriimmten Teil 
fiir kleine Variationen die Neigung der Tangenten in den Punkten D, bzw. D, 
an die Kennlinien K, und K, massgebend. Wiinscht man in dem sogenannten 
Nennbetriebsgebiet, also bis zum Nennstrom J, und bis zur Nennspannung U,, 
tiberall parallele Gerade als Kennlinien zu erhalten, so lasst sich der hiefiir not- 
wendige gerade Teil von K, bzw. K, sofort abgrenzen. 

Das dynamische Verhalten kann an Hand der Blockschemas Figur 4 und 5 
beurteilt werden. Zu der unabhangigen Variablen w,(¢) tritt der willkiirliche 
Belastungsstrom 7(¢) als «Stérgrésse». Der Momentanwert u(¢) der Klemmen- 
spannung kann zum Beispiel als Uberlagerung einer Exponentialschwingung 
mit der komplexen Frequenz p =o + 7 @ (j2= —1) iiber einen konstanten 
Wert U aufgefasst werden: 

u(t)= U + Re(U e??) . (3} 


Im linearen Teil der Kennlinien K, und Ky gilt dann nach Figur 5 fiir den 


Zeiger U: 


U = 6,(6) 0, = Gey 7, (4) 
2k eee bee 
OAV Meader aote +Ke(p) 14+ T/1+ K g(p)]p’ 


(4a) 
bo See ee oe 
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Las LP lp S 
rrr tap ee ee | 
: : : (4b) 
os agrees errr | 


_ R, beziehungsweise Z, ist das die Spannungsregelung bedingende Glied. Obige 


_ Formeln fithren fiir = 0 und g(0) = 1 zum bereits behandelten statischen Ver- 
“halten nach Gleichung (2). 


= 


7 


Figur 4 
Blockschema des geregelten Generators. 


Figur 5 


Resultierendes Blockschema. 


Das dynamische Verhalten kann fiir den Fall g(f) = 1, wenn also der Ver- 
starker viel rascher arbeitet als die Gleichstrommaschine, gut tibersehen wer- 
den. Fiir Variationen U, der Bezugsspannung geht die relativ grosse Zeitkon- 
stante 7, des Feldes nach (4a) auf einen (1 + K)-mal kleineren Wert zuriick. 
Dasselbe tritt fiir die Quellenimpedanz Z, bei Stromvariationen J ein, jedoch 
nur fiir ganz kleine Werte von #; fiir grosse nahert sich der Bruch rechts im 
Ausdruck fiir G,(p) der Einheit. Die geregelte Maschine reagiert also rascher 
auf Verdnderungen als die ungeregelte. Eine plotzliche Belastungsschwankung 
kann durch eine sprungweise Variation der komplexen Impedanz Z erfasst 
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werden. Mit der reziproken komplexen Ubertragungsfunktion 


h,(p) = 76), (4c) 
des Verstarkers wird F $ 
1 
Cl?) = TER) ALE’ es | 
Z hy(p) 2 
Gold) = KT HO) IE a 


Der Ausdruck fiir G, zeigt die Reduktion der Quellenimpedanz fiir p = 0, der 
Term rechts gibt allgemein den Frequenzgang von =0 (K:1+ K) bis p=oo(K). 

Figur 6 erlautert die Rtickwirkung der Netzimpedanz Z auf die Spannungs- 
haltung und fiihrt zu 


ee. 3 —— 5 ret — ae . Af 
k 1+ h,() (@/K) [1 + (2,/2Z)] me See 


Die technisch wichtigsten Spannungsschwankungen treten bei konstantem Uy 
infolge von beliebigen zeitlichen Variationen von Z auf. Leider kénnen die 
dadurch entstehenden nichtlinearen Differentialgleichungen durch die Laplace- 
Transformation nicht erfasst werden. 


Figur 6 


Blockschema mit Belastungskreis. 


Fur g(p) = 1 kann der Verstairkungsgrad K beliebig hoch gewahlt werden, 
ohne dass Instabilitat auftritt. Hat jedoch der Verstarker im Vergleich zur 
Maschine eine wesentliche Zeitverzégerung, wie solche zum Beispiel bei magne- 
tischen Verstarkern auftreten, so miissen zur Beurteilung des maximal zulas- 
sigen Wertes von K die Stabilitatskriterien herangezogen werden. 

Regler fiir konstanten Strom lassen sich analog graphisch und dynamisch 
behandeln. Die Quellenimpedanz wird dabei (1 + K)-mal grésser, die dusseren 
Kennlinien des geschlossenen Kreises fallen steil ab. Analog lassen sich Kaska- 
denschaltungen von mehreren Maschinen oder Verstarkern und von vielen 


Arten von Proportionalreglern unter Beniitzung ihrer Kennlinien graphisch 
behandeln. 


i 
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B. Antrieb fiir konstante Drehzahl 


; Die Kontrolle der Drehzahl des Motors M und der Arbeitsmaschine A 
4 geschieht nach Figur 7 durch einen gesteuerten Stromrichter und mittels des 


Figur 7 
Gesteuerter Stromrichterantrieb mit Tachometerdynamo und Verstarkerroéhre. 


Tachometergenerators 7, dessen Gleichspannung U; mit dem einstellbaren 
Sollwert U, verglichen wird. Die folgenden Uberlegungen gelten analog, 
wenn statt des Stromrichters ein Gleichstromgenerator verwendet wird. Zu- 
nachst sollen die statischen Verhaltnisse ahnlich wie im ersten Beispiel graphisch 
behandelt werden. Der Fehler (E = U, — Ur) wirkt auf das Gitter der Ver- 
starkerréhre V, deren Ausgangsgleichspannung U, dem Gittersteuergerat St 
zugefiihrt wird. Die Veranderung der Leerlaufgleichspannung U,, am gesteuer- 
ten Stromrichter wird in bekannter Weise durch zeitliche Verzégerung des 
Ziindeinsatzes der einzelnen Anoden gegeniiber dem nichtgesteuerten Zustand 
erreicht. Beim m-Phasen-Gleichrichter, zum Beispiel m = 6, sinkt dann die 
Leerlaufspannung U, gegeniiber ihrem Maximalwert U,, ab laut 

(pee COSé (5) 


lo 


mit « als Ziindverzégerungswinkel. 
Die dusseren Kennlinien des gesteuerten Stromrichters (Figur 8) sind paral- 


lele Gerade ahnlich denjenigen einer fremderregten Gleichstrommaschine (Fi- 
gur 2). Das Steuergerat S¢ erzeugt nun einen um so grésseren Steuerwinkel «, 
je negativer die Eingangsspannung U, ist. Die Kennlinie K; in Figur 8 stellt U, 
_in Funktion von U, dar. Aus dem Kreis K mit dem Radius U,, lasst sich fiir 


ZAMP V/29 
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jeden Wert U, der Steuerwinkel « sowohl fiir Gleichrichter- (0 << « < 90°) wie © 


auch fiir Wechselrichterbetrieb (90° < « < 180°) konstruieren. 


hic cHiees 


Der Stromrichter weist einen ohmschen (R,) und einen induktiven, durch 


die Wirkung der Streuung L, wahrend den Uberlappungszeiten bedingten 
Spannungsabfall auf. Unter Beriicksichtigung des Ankerwiderstandes R, er- 


Figur 8 


Statische Kennlinien des Stromrichterantriebes. 
geben sich daher fiir den Mittelwert U,; der induzierten Ankerspannung laut 


U; =U, - (Ry + Rat 4% fL,)I=U,— RI (6) 
parallele Gerade (Figur 8, AB|| LN) als Kennlinien (f = Frequenz). Bei vielen 
Problemen der Servotechnik erweist sich die Einfithrung von bezogenen Grés- 
sen, also von homogenen Koordinaten, als sehr zweckmissig. Als Bezugsgréssen 
werden die «Nennwerte» gewahlt, also der Nennwert der Drehzahl n,, des 
Stromes J,,, des Drehmomentes M,, des Motors und die maximale Gleichspan- 
nung U,,. Die bezogenen Werte lauten: 


mM M, : o 
ate M,, ee = M,’ age CR Sean pe (7) 


n nr m 


wobei M, das Gegendrehmoment der Arbeitsmaschine bedeutet. Fiir die vom 
konstanten Feldfluss ® induzierte Ankerspannung U, und fiir das Drehmoment 
M gelten bei Vernachlassigung der Ankerriickwirkung: 


U;=C ow, U,= Coa, , M aCigiMe= Cie (8) 
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und daher ist: 


U; n wo M 1 
y= ~ == = — 
on Ny Me™ or M,, iE (8a) 


Durch Division der Gleichung (6) mit U,, folgt der Zusammenhang zwischen 
der bezogenen Drehzahl » und dem bezogenen Drehmoment yu im offenen 
Regelkreis (Figur 9) unter Vernachlassigung der Reibung zu: 


Tig 
Vi We i A ee (9) 

m 

7 

dy, 
ie 
, Figur 9 Figur 10 
Bezogene Drehzahl-Drehmomenten-Kennlinie Bezogene Drehzahl-Drehmomenten-Kennlinie 
des Stromrichterantriebes. der Arbeitsmaschine. 


Die entsprechenden Kennlinien der Arbeitsmaschine zeigt Figur 10 fiir die 
Falle: 
I: “a = konstant (zum Beispiel Hubwerk), 
Tes, eg, (near), 
III: uw, = v2 (zum Beispiel Pumpe, Ventilator), 
IV: beliebig. 


Beim offenen Regelkreis fallt nach Figur 8 die induzierte Spannung U; und 
damit die Drehzahl » nach der Geraden AB ab, falls AyB = RI gewahlt wird. 
Wie liegen nun die Verhdltnisse beim geschlossenen Regelkreis? Nimmt man 
fiir einen vorgegebenen Strom J zunachst einen wiliktirlichen Spannungsabfall 
AC, an, so ergibt eine zu Figur 2 analoge Konstruktion mittels der Kennlinien 
K, und Ky des Verstarkers beziehungsweise des Steuergerates eine Erhohung 
der Quellenspannung U, um A C,. Die Verbindungslinien C,C, D,D usw. fir 
verschiedene Werte von AU gehen nun alle im linearen Bereich durch einen 
Punkt S. Es stellt sich diejenige Gerade D,D im geschlossenen Regelkreis ein, 
welche parallel zu den Kennlinien fiir U;, also parallel zu AB, liegt, woraus 
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unter Vernachlassigung der Reibung folgt: 
AU, = KAU; AU AU (1+ K)AU=RI1, 


Ur eke ath in 

Cass Le Re K, Pi i ( ) 

Uae ly eek eee (11) 
WARK 2 mg $s 


Beim geschlossenen Regelkreis ist also der Drehzahlabfall (1 + K)-mal kleiner 
als beim offenen. K ist der Verstarkungsfaktor im geschlossenen Regelkreis. 


Fiir die Untersuchung des dynamischen Verhaltens schneiden wir den Re- 


gelkreis an den Steuergittern des Stromrichters auf (Figur 11), fiigen dort neben 


Figur 11 
Aufgeschnittener Regelkreis. 


einer konstanten Spannung eine St6rspannung ein und untersuchen das Ver- 
halten aller Regelgrdéssen unter der Wirkung dieser erzwungenen Schwingung, 
die wir als eine allgemeine Exponentialschwingung von der Form 


Re [A e?']= Ae” cos(wt+ q) 


mit der beliebigen Kreisfrequenz w, der beliebigen negativen oder positiven 
Dampfung o, also von der beliebigen komplexen Frequenz =o +7 w wahlen. 
Die massgebenden Gleichungen folgen einerseits aus dem Spannungsgleich- 
gewicht des elektrischen Kreises und anderseits aus dem mechanischen Kreis 
nach NEWTON. 

Die am inneren Widerstand R, und an der Streuinduktivitaét L, des Strom- 
richters, an einer eventuellen Glattungsdrosselspule Ly, am Ankerwiderstand 
R, und an der Ankerinduktivitat L, auftretenden Spannungen bedingen eine 
Spannungsdifferenz zwischen den Momentanwerten u,(¢) der Quellenspannung 


und der induzierten Ankerspannung w,(), welche (angenahert infolge Anderung 
der Uberlappungszeit) betragt: 


20% : } : ; 
u(t) — u,(t) = (R, +R + wre ae | tae tds aerate Tan) - = ee soo S : 
(12) 
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_ Fair die komplexen Zeiger folgt daraus: 
4 


U, — O,= [(Ry + Ra+ 2% fL,) + (Ep + Lat L,) A] T, | 


tg Rot Rot =" 1, L=1,+1a+lLy,, (13) 
U,-U=(R+pL)T=RA+7AT, B= = | 
und daraus durch Division mit U,, fiir die bezogenen Zeiger: 


—~ —_ RIn i s RL 
Ge dai om (l+ ZA) z =v (l+ Ta) ph, y= fo (14) 


_ Die Differenz der Drehmomente M(t) des Motors und M,,(¢) der Arbeitsmaschine 
_ wird fir die Beschleunigung der Massen und fiir die zur Drehzahl proportional 
_ angenommene Reibung aufgebracht laut 


Mi) —M,) = J 2247 M,, (15) 


woraus durch Division mit M,, fiir die komplexen Zeiger folgt: 


Be = Tee bI+ 1D = vv +T, pr, (16a) 
ey aT fet wenn, (16b) 
je — Ha = (+ Td) D=7(L 4 Im) P= 7 Ep? (16¢) 


T,, ist die Zeit zum Hochfahren der Maschinen unter der Wirkung des Nenn- 
_ drehmoments (M, = 0) bis zur Nenndrehzahl ,. Linearisiert man nach Figur 10 
fiir kleine Variationen den Zusammenhang zwischen Drehzahl und Drehmo- 
ment der Arbeitsmaschine laut 


Ha = Mag t+ Ka OW, Ky= a Ha= Kara, (17) 


a 


_ so folgt schliesslich fiir die starre Welle (v, = v) beim offenen Regelkreis: 


far (l+T,.p)> mit r=r+Ki, In=—a Tn- (18) 


y 


Die in diesen Gleichungen auftretenden Gréssen, wie die Impedanz (Kk + p L), 
die «bezogenen» Impedanzen oder Zeitglieder 


2,2 L472; t= 1S Ts. Pp (19) 
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des elektrischen bzw. mechanischen Kreises sind in den Figuren 12a, bund c — 
fiir komplexe Werte von p = | p | exp (j y) ¢ dargestellt. Aus diesen ergeben sich 
die ohmsche und induktive Spannung RI bzw. pLT, das Reibungsmoment — 
M,,7 vin Phase mit » sowie das Newtonsche Beschleunigungsmoment M,T,pyv~ 
um y voreilend. Die Differenzspannung (U,, — U;) bewirkt den Ankerstrom ve | 
dieser das Drehmoment M; dieses, vermindert um das Gegenmoment M,, be- — 


opens 
\ 
> 7 
De J 


Z-/+yrZe2m 


Zz Wy 
Yen 
: = \ Gane (B+ Gr) 
0 U; M,rv Tang 1 or 


Figur 12 


Zeigerdiagramme ftir eine komplexe Frequenz p. 


wirkt die Drehzahl x, die ihrerseits die induzierte Spannung U, hervorruft. 
Nimmt man U, oder ¥, als Ursache, » als Wirkung, so findet man aus der 
Figur 12c die Zusammenhinge 

U,—U,=R2z,1, M—M,=7M, 2m? (20) 
und damit die fiir die kombinierte Wirkung des elektrischen und mechanischen 
Kreises charakteristische inverse Ubertragungsfunktion 2(p) zwischen den be- 
zogenen Leerlauf- und Lastdrehzahlen 
30) = = 1 ty 1,3, (21) 


v 


&. 
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_ Man sieht daraus, dass man diese Zusammenhiange graphisch auch fiir komplexe 
_ Frequenzen darstellen kann, wahrend dies bis anhin nur fiir ungedampfte 
_ Sinusschwingungen (p = 7 w) tiblich war. z() vermittelt als quadratische Funk- 
_ tion in # eine konforme Abbildung der «Frequenzebene p» auf eine zweiblattrige 
_ Riemannsche Flache 2(p). Das Bild der imagindren Achse (p = 7 w) ist eine 
_ Parabel und wird als komplexe Frequenzgangkurve bezeichnet (Figur 13) 


47 p-Ebene 


Figur 13 Figur 14 


Zeigerfrequenzgangkurve der kombinierten Blockschema ftir die bezogene Drehzahl v. 
Impedanz z. 


_ Eine Exponentialschwingung » bewirkt tiber das Tachometer Figur 11, itber 
die Verstarkerréhre V und das Steuergerat St eine Exponentialschwingung »¥* 
- der Leerlaufspannung des Stromrichters im aufgeschnittenen Regelkreis. Es sei 


: - ot K 
i = =i *s. e5(6) 8,(6) ¥ = —K 80) ¥ = — pa? a) 


mit g(p) als resultierendem Zeitglied. Eigenschwingungen des geschlossenen 
' Regelkreises resultieren aus »¥ = ¥,, es sind dies ebenfalls Exponentialschwin- 
gungen, deren komplexe Frequenzen daher die Bedingung 


z(p) + K g(p) =0, 4,(p) 2(p) + K=0 © (23) 


erfiillen. Das System ist stabil, wenn alle Wurzeln dieser Gleichung negative 


Realteile aufweisen. - . 
Die bezogene Drehzahl y lasst sich im linearen Bereich als Funktion der 


zwei unabhangigen Variablen U, und M, nach Figur 14 darstellen: 


¥ = Gy(p) Uy — Gal) Ma (24) 
mit 
K e(P) x gee 24 
GP) =z ap) + Kee) Kr wMpz+ kK’ see 
Zs ‘ h(P) Ze 24b 
CP) =1 sp) + Ke)’ WP)Z+K™ ae! 


5 

i 
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G,(p) erhalt man durch mehrmaliges Anwenden der Formel : 


fiir den gegengekoppelten Kreis nach Figur 15 auf Figur 16 mit yw, als Ursache | 
und bei konstantem Up. Dieses Verfahren wird bei komplizierten Servosystemen ~ 
mit mehrfachen Gegenkopplungen mit Erfolg angewandt und entspricht dem 
Auflésen von Gleichungen mit mehreren Unbekannten. ; 


Ms 3 v 
UP Pee 


Figur 15 Figur 16 


Gegengekoppeltes System. Blockschema zur Ermittlung der Ubertragungsfunktion 
zwischen M,(t) und q(t). 


Sind die Zeitverzdgerungen des Verstarkers und des Steuersatzes vernach- 
lassigbar gegentiber den tibrigen, ist also g(#) = 1, so gehen die Eigenfrequenzen 
aus der Gleichung 


z(p) + K =0 (26) 


hervor. Hiefiir fallt in Figur 13 P auf Ky, und es sind zum Bildpunkt K, die 
zwei Originalpunkte in der #-Ebene zu suchen. Denkt man sich den Verstir- 
kungsfaktor K von Null ansteigend, so interessiert das Wandern der zugeh6- 
rigen Wurzeln ,, p, in der Frequenzebene. Diese Wurzelmethode kann zur 
Bestimmung der oberen Grenze von K bestimmt werden, die eintritt, wenn 
Wurzeln die imaginare Achse iiberschreiten, das System also instabil wird. Man 
sucht dazu das Bild der negativen reellen Achse der z-Ebene in der -Ebene 


Sd 


% 
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auf. Im vorliegenden Falle liegen die Verhaltnisse besonders einfach. Aus (21) 


folgt 
aq 12 7G 


4 fein PS Te Te) be peel (27) 
und daraus 

1 1 1 : 

amines (Gee 7) (28) 


also liegen komplexe Wurzeln unabhingig von K stets auf der Senkrechten im 
‘Abstand (— oc) vom Ursprung. Reelle Wurzeln sind zudem negativ. Sind die 
beiden Wurzeln des offenen Regelkreises beispielsweise negativ reell, und ver- 
starkt man nun die Gegenkopplung sukzessive durch Steigern von K, so gehen 
die beiden Wurzeln zunachst in eine Doppelwurzel (—o) iiber und werden 
daraufhin komplex mit zunehmender Frequenz w (Figur 17). Je starker die 
-Gegenkopplung, um so hoher ist also die Frequenz der Exponentialschwingung 


ef 
, 


O —=> periodisch 
in 


Figur 17 Figur 18 
Das Wandern der Wurzeln in der komplexen Die drei Gebiete fiir aperiodische und 
E Frequenzebene. periodisch gedampfte Ubergange. 


bei gleichbleibender Dampfung o. Ein stabiles System zweiten Grades wird 
infolgedessen durch eine beliebig starke, keine Zeitverzégerung aufweisende 
Gegenkopplung nicht instabil. Dieses Wandern der Wurzeln bei wachsender 
Gegenkopplung kann analog studiert werden fiir Systeme hoheren Grades in #. 
Man geht vom offenen Kreis aus, wobei fiir stabile Systeme alle Wurzeln auf 
der negativ reellen Achse liegen und bei zunehmendem Wert von K eventuell 
tiber die imaginare Achse in die rechte Halbebene tibertreten kénnen. 

Die wichtige Grenze zwischen aperiodischem und periodisch gedimpftem 
Vorgang, also bei der kritischen Dampfung mit verschwindender Determinante, 
fiihrt bei g(p) = 1 zu den zwei Gleichungen 

1 in 4 Wh 2 
e = z (1 ; yrt +] aoa ET. oe (l+pyr+ | (1+ A). 
“ (29) 


Diese Beziehung zwischen den beiden Zeitkonstanten T,, T,, wird durch zwei 
Gerade nach Figur. 18 dargestellt. Schwingungen zwischen dem magnetischen 
Energiespeicher L und dem mechanischen Energiespeicher J sind um so schwie- 
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riger zu vermeiden, je “starker die Gegenkopplung ist. Um so kleiner wird 
gleichzeitig die Periodendauer der gedampften Schwingung und erreicht schliess- 
lich einen Wert, fiir den die Voraussetzungen der Theorie nicht mehr giiltig 
sind. Denn beim 6-Phasen-Stromrichter bei 50 Hz Netzfrequenz erfolgt die 
Wirkung der Steuergitter unstetig nach je 3,33 ms, die Signale weisen also 
eine «Totzeit» auf. Zu hohe Verstarkungen sind also zwecklos. Glattungs- 
drosseln L, erhéhen T, und damit die Neigung zu Schwingungen. ; 

Dem folgenden Zahlenbeispiel liegt ein plétzliches Abfallen des Gegendreh- 
momentes um den Wert AM,, laut m,(t) = + Af1}, zugrunde. 

Die Laplace-Transformierte von y fiir diese Schrittfunktion (als Spezialfall 
einer Exponentialschwingung # = 0) von ju, ergibt sich bei g(p) = 1 aus (24b) 


a Pe 
Liv) =yA—-= oat 30 
”) p26) + K “4 
deren Lésung fiir reelle Wurzeln #, f, lautet: 
yp — 
erg eer a tg (31a) 
und fiir komplexe (y’ = a — y): 
us we — oe —ot ‘sin (@ t+ y’) 1 Sot ° 
LOT yy (1 e a )+ i ‘sinwt. (31b) 


Die Figur 19 zeigt den Verlauf von » fiir den offenen (K = 0) und den geschlos- 
senen Regelkreis mit K = 10, wobei fiir letzteren der Endwert 11mal kleiner ist 


I, =1s, T,=10ms, y = 0,1, r= 0,02, T,, = 50s, o = 50,01 Hz. 


m 


@m = 92,1 Hz, f= 14,6 Hz, T= 68,3 ms. 


Figur 19 


Oszillogramm der bezogenen Drehzahl y(t) beim offenen und geschlossenen Regelkreis. Lage der 
Wurzeln in der Frequenzebene. 
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und zudem viel schneller erreicht wird, da die Anfangstangenten fiir beide 
Falle identisch sind. Darin liegt der Grund fiir die « Verkleinerung der Zeitkon- 
‘stanten» bei steigendem K. Nimmt man als Zeitkonstante die Zeit bis zur Er- 
reichung des Endwertes mit der Anfangsgeschwindigkeit, so verhalten sich 
‘diese wie 1:11. 

Haben der Verstarker und der Gittersteuersatz wesentliche Zeitverzoge- 
rungen [g(f) + 1], so kann das System instabil werden. 


i C. Drehzahlregelung eines Walzmotors 


Das folgende Zahlenbeispiel!) behandelt nach Figur 20 einen mittels einer 
Tachometerdynamo 7, eines Verstaérkers und einer Amplidynmaschine gere- 


bersl8rker Amp lidyn Generator Rt, Motor Arbeitsmaschine 
heme he } M(t) M(t) 
Vu 


Tachometer- 
dynamo 


Figur 20 
Schema des Walzmotorantriebes. 


gelten Ward-Leonard-Satz zum Konstanthalten der Drehzahl eines Walzmo- 
tors von 


Boa 183 kW Nennleistung, 
tne Z5U A Nenngleichstrom, 
U,,= 628 V induzierte Nennankerspannung, 


Ny = 150 U./min = 2,5 Hz Nenndrehzahl. 


Wir betrachten sukzessive die Ubertragungsfunktionen der einzelnen Organe 


Figur 21a). 
a Verstarker 


Ky = 12,2). = 1. 


Amplidynmaschine : 
Langsfeld: Ry= 802, 1, = 8H, pul s,* 4 = 14+ 7,2; 
Querfeld: R,= 1,062, L,=0,318H, 7,=0,38, %=1+ T, p; 
Maschinenkonstante: C, C, = 287022, voll kompensiert. 


1) Harotp CuestNut und Ropert W. MAYER, Servomechanisms and Regulating System 
Design, Bd. I (John Wiley and Sons Inc., New York und London 1951), Seite 468. 
med, 
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Generator : 

R, = 7,152, L,=— 2440, <= 315s Coa 1B e Za =14T,; p. ; 
Motor 

Induzierte Spannung: U;=Cw, U;,=Ca,, 
Drehmoment : ee Om a: | Rea es beg 
Maschinenkonstante: C =40NmA-1=40Vs, : 
Winkelgeschwindigkeit: w, = 22, = 15,7 Hz.. H 


Stromkreis — Generator — Motor 


Widerstand: R,=12m2, 
Induktivitat : L[,=0  vernachlassigt. 
Arbettsmaschine 


M,=konstant (Fall I von Figur 10), g,(p) = 0. 


Aggregat Motor-Arbeitsmaschine 


Tragheitsmoment: J = 6000 kg m?, 
Reibung: ge 


Tachometer 
Us a Kr Y, Kr == 15,708 ‘Ve 


In Figur 21a sind die Ubertragungsfunktionen der einzelnen Glieder darge- 
stellt. Das Zeitglied des aus beiden Ankern bestehenden Kreises ist 


a=l+T,o mit Dips LAs 
Aus Figur 216 erhalt man: 


v = Gy(P) Up — G2(P) fas (32) 
GO) > ee apy eer OO =7 Saas eae ee 
mit 
221 pp Py = 2. (32b) 
BP) = 8 Fe Bon (32c) 
K a om z z = 24,236. (324) 


Mit der reziproken Ubertragungsfunktion h,(p) des Verstarkers 


hy(p) = ae (32e) 
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erhalt man: 
Gi) = 3 — Gama 33 
Kr | h,(6) hip) + K’ ae 
: Gy = y 22. hel) 0) 
oR Gh) ap) RE ul 
ann Amplidyn . Generetor Motor Arbeitsmaschine 


slarker: 


Figur 21 
Blockschema des Antriebes nach Figur 20. 


Im vorliegenden Beispiel wurden L, und y vernachlassigt, damit erhalt man: 
go), ol ey t= mL pH lt p= %, 
Oy J R t 
Ig=Yy Lena Ae = - aw = 45 ms (34) 
und fiir die Hochfahrzeit T,, = 1,88 s. Der Ausdruck z ist also in eine lineare 
Funktion von # tibergegangen. Der Spannungsverstarkungsfaktor zwischen der 
Fehlerspannung F und der Quellenspannung U,, des Generators betragt 
Uy, _ Ky Oy C2 Ca 


Ks se 


: Bop ess IOI AG. (35) 


Die Ubertragungsfunktionen G, und G, fiir den offenen Regelkreis sind beson- 
ders einfach : 


K, 1 
G,(p) = U,; hyp) h(p) ' ¥ 
Za 1 


Gl) =y Baye. (37) 
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Slide 


Das Produkt h(p) der vier Zeitglieder ist in der Kenneleyschen Schreibweis 
h(p) = Bo 21 22% = (1+ Ty) (1+ Fp) (1+ 726) (1 + 158), 
BO) = ET Ag 2 Pe cage %,=14+ 7, p=A4, /9,- 


ee 


Die Eigenfrequenzen des gesamten Systems ergeben sich aus den Nullstellen von 


hy(p) h(p) +K = 0. . " (38), 


Der fiir viele Servoprobleme charakteristische Ausdruck h,(p) h(p) ist oft ein. 
Produkt von linearen und quadratischen Gliedern in #. Ist speziell g(p) = 1, 
so wird ; 


h(p) = —K. (39) 


Wenn die vier Wurzeln dieser Gleichung negative Realteile haben, so ist das 
System stabil. Will man die Bestimmung dieser Wurzeln vermeiden, so muss 
man eines der Stabilitatskriterien anwenden. Wir zeichnen beispielsweise das 
Bild der Geraden = 7, also die komplexe Frequenzgangskurve h(j w) in 
der h-Ebene, Figur 22a (Leonhard-Kurve). Die rechte #-Halbebene (Figur 22d) 
geht dabei in das schraffierte Gebiet (Figur 22c) tiber. Da nun der Punkt 


p-Lbene h-Ebene 


a 


Figur 22 
Leonhard-Kurve der Funktion h(p). 


h(j w) = (1 + 0,045 7 w) (1 + 0,17 w) (1 + 0,37 w) (1 + 3,157°w) = Hel® 


)( 
a) 0,3 1 | 2 3 + | 5 | 10 | 
A | 1,38 | 3,46 | 7,60 | 13,4 | 21,3 | 33,6 15452 
7) 51°04’ 97° 24’ 128° 26’ | 150" 22° || 16724344) 181°53¢ 229° 


Soy 


ey 
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(—K = —24,24) im nichtschraffierten Gebiet liegt, ist das System stabil. Der 
hochstzulassige Wert des Verstarkungsfaktors K der Schleife liegt bei 33,2, 
wenn der Punkt K auf den Schnittpunkt S fallt. Die Marge ist also nicht sehr 
‘gross. Will man vorsichtig sein, so wird man ein stabilisierendes, phasendre- 
hendes Glied einbauen. 

Die Nyquist-Methode verlangt nach Gleichung (33a), (330) oder (38) das 
Aufzeichnen der Kurve [K: h,(j w) h(j w)], wofiir in den USA Vordrucke ver- 


HTP 
‘ LAS 
a8 wee, le XO 


b 


Figur 23 
Nyquist-Diagramm. 


wendet werden?). Befindet sich wie in den Figuren 23 und 236 der Punkt 
(— 1/07) ausserhalb der schraffierten Flache, so liegt Stabilitat vor. Figur 24 
zeigt schliesslich die doppeltlogarithmische Darstellung der Amplitude H und 
Phase y von A(j w) sowie den zugehGrigen charakteristischen Polygonzug der 
Asymptoten nach Bone. Man sucht hier den Schnittpunkt S fiir yg =z und 
ermittelt, ob der zugehorige Amplitudenwert H(S) grésser als K ist. Oder man 
sieht zu, ahnlich wie in Figur 22, ob die Phase m; des Punktes T fir H = K 
kleiner als z ist. Le 

Noch einige Zahlenwerte zur Erlauterung der Prdzision. Die Quellenspan- 
nung des Generators betragt bei Nennbetrieb R,/, = 15 4 mehr als die indu- 
zierte Motorspannung U,; = 628 V. Im offenen Regelkreis steigt deshalb bei 


1) macoun Cuestnut und Ropert W. Maver, Servomechanisms and Regulating System 
Design, Bd. I (John Wiley and Sons Inc., New York und London 1951), Seite 468. 
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voller Entlastung die Drehzahl um 


Die Spannung U, = E ist dabei K, = 969,46mal kleiner und betragt 663,3 mi 
Wie sich aus (32) und (33b) fiir du, = 1 und U, = 0 entnehmen asst, betragt 
im geschlossenen Regelkreis bei voller Entlastung die Drehzahlerhohung . 


SB ibis nie gO al z 
a a ae ri 230%: 


nt 


fet aes ib liake 
Avs tie =O, Dai 


rye iat pa nab ty ea ali 


mens 


| 


(40) 


wir haben also eine Ein-Promille-Steuerung. Bei voller Entlastung von der 
Nennleistung steigt die induzierte Spannung von 628V um 0,95°/,, auf 


o 


Z z 
2 G GL 
0 


0 


=a0 


Aju) =H./p 


SOL 
OF O4 O5 


° 


Figur 24 


Bode-Diagramm. 
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628,594 V, wahrend gleichzeitig die Quellenspannung U,, von 643 V auf diesen 
Wert zuriickgeht. Die Tachometerspannung U; steigt ebenfalls um 0),95°/o9, Vou 
15,708 V auf 15,723 V. Diese minime Anderung um nur 14,9 mV bewirkt also 
die totale Anderung der mechanischen Leistung um 785 kW. Dabei bleibt die 
Vergleichsspannung konstant bei U) = 16,371 V, wahrend die Fehlerspannung 
von E,, = 663,3 mV auf Ey = 648,4 mV zuriickgeht. Sie ist ja proportional zu 
Uys, aber 969,46mal kleiner. Man kann daraus ermessen, welche hohen Anforde- 
Tungen an die Prazision und Konstanz aller Teile des gesamten Regelkreises 
gestellt werden. 

_ Ist die Drehzahl vy = y, nach den Figuren 10 und 216 durch K, £qa(p) mit dem 
Drehmoment yz, verkniipft, so ergibt sich aus (32) eine zu (4c) analoge Bezie- 
hung 

G,(P) 


* T+ [ye Ka ea(b)e + K elP)] * 


Andert sich K, g,(p) zeitlich, so resultiert eine nichtlineare Differentialglei- 
chung. 

Die experimentelle Untersuchung des Frequenzganges von Servosystemen 
umfasst gegentiber der Rechnung weitere Einfltisse, wie Wirbelstréme, die ma- 
gnetische Hysteresis, die Reibung usw. Wie man aus den Werten von w der 
Figur 24 ersieht, bendtigt man hiefiir zur Untersuchung grésserer Aggregate 
Generatoren mit sinusférmiger Spannung und enorm kleiner Frequenz, zum 
Beispiel 0,1... 10 Hz. 

Soll die Eisensattigung mit beriicksichtigt werden, so miissen zur Behand- 
lung der dadurch auftretenden nichtlinearen Differentialgleichungen Analogie- 
rechengerate mit nichtlinearen Gliedern herausgezogen werden. 


Résumé 


Cette étude traite trois exemples tirés de l’électro-servotechnique. On construit 
d’abord dans un graphique les caractéristiques extérieures d’un générateur réglé 
en circuit fermé. Puis on discute le réglage automatique du nombre de tours 
d’un moteur commandé par un redresseur. Enfin le réglage en circuit fermé d’un 
moteur accouplé a un laminoire est traité par un calcul numérique ; les différents 
critéres de stabilité permettent de déterminer l’amplification de limite. 


(Fingegangen: 31. Marz 1954.) 
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Deformations Possible in Every Isotropic, Incompressible, 
Perfectly Elastic Body 


By JERALD LAVERNE ERICKSEN, Washington, D.C., U.S.A.) | 


a meee 


1. Introduction ; 


The problem considered here is that of determining deformations which can 
be produced in every isotropic, incompressible, perfectly elastic body by the 
application of surface tractions when there is no body force acting. Examples” 
of such deformations may be found in the general solutions for homogeneous 
deformation [1]?), [2], torsion of a right circular cylinder [1], [3], bending of a 
block [4], deformation of a spherical shell subject to internal and external 
pressures [5], and in a five-parameter family of solutions recently obtained for 
materials possessing transverse isotropy [6]. These latter solutions were disco- 
vered in an attempt to solve the problem discussed here. Three of the solutions 
obtained here are new. 

Let-the points of an undeformed elastic body be assigned coordinates X* 
in a coordinate system with metric tensor’) G, ,. Suppose the body is subjected 
to a deformation, and let the points of the deformed body be assigned coordi- 
nates x* in any independently selected coordinate system with metric tensor g; ;. 
Greek indices always refer to the coordinates X*, Latin to the coordinates x’. 
It is assumed that the deformation is described by a one-to-one mapping, 


IF =a LX (1.1) 


which is at least three times continuously differentiable. 
When the inertial and body forces vanish, the equations governing the 
behavior of a material of the type considered here become 


a (1.2) 
: 4 Pay _1W\t 0x j 
We RO 8 ae Oe tee ts 
fers di, (1.4) 


) Applied Mathematics Branch, Mechanics Division, Naval Research Laboratory. 
) Numbers in brackets refer to the Bibliography on page 489. 


1 
2 
3) We use standard tensorial notation. See, e.g., [7]: 
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where ¢ is the stress tensor, # is an arbitrary pressure, 
_ of the matrix || c;|| of the tensor ci, 


ci| is the determinant 


tent Oe Oe 
C5 Cup Oxi . Axi? (1.5) 
-1\ij _ fap Oxt Oxt | 
(c ) a G OxXxe 4 OXB’ (1.6) 
Sve 1 SUN Ve 
L=(, Was [P— evi e-4,], (1.7) 


and & = 2'(I, II) is the strain energy, assumed to be at least three times con- 
tinuously differentiable with respect to I and II. As is easily verified, 


(aay? Ca i . (1.8) 


Since we are working in Euclidean space, the Riemann tensor based on G, , 
vanishes. It then follows from (1.5) that the Riemann tensor R;;,, based on 
c;; must vanish. Analytically?), 


or <2 ( ee Oi 
ape? Ox) oxk Ox? ox! Ox) Ox! Ox* oxk 
0c; OCp OC; 0c; Oc 0c; 
ae co mn ( ym , uae +i i] ( mn a 25 a tl 
og Oxk Ad Oxi on™ ox? a Ox Oxt 
(1.9) 
Sa ( OCj m sf Om OCj1 
Ox! Ox) oxnm 
OC; n OCk n OC; x 
x (— + = =0. 
( Oxk 3 Ox Ox” ) 0 


For R;;,; to be the zero tensor, it is necessary and sufficient that the six com- 
ponents Ryo12, Reso, Rsisi, Rives, Rossi, Rsr1z Vanish. 

Suppose we are given a symmetric tensor c;; such that the matrix || c;;|| 
is positive definite. We can then use (1.8) to determine (c~1)'7, If these tensors 
satisfy (1.9), and if G,, is a given Euclidean metric tensor, there will exist, at 
least locally, a mapping of the form (1.1) such that (1.5) and (1.6) are satisfied. 
In order that c,; determine deformations?) possible in an incompressible mate- 
rial, it must satisfy (1.4). Suppose these kinematical conditions are satisfied 
and that the function (I, II) is specified. Then, if it is possible to determine ¢ 
and p:suchsthat (1.2) and (1.3) are satisfied, c;; corresponds to a deformation 
possible in a material for which & is the strain energy. Even if such a deter- 


1) TRUESDELL has compiled a long list of papers in which equations (1.9) or equations equiv- 


alent to these are discussed. See Math. Rev. 12, 556 (1951). 


2) As is well known,:c,;; determines the deformation only to within a rigid motion of the un- 


-deformed material. 
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mination be possible for one choice of X, it will not, in general, be possible for 
arbitrary choices of Y. We will determine the exceptional tensors ¢,;; such that 
this determination is possible for every choice of 2 under rather weak restric- 
tions on the form of c;;. We shall compute the corresponding stresses and — 
describe a simple physical situation to which the solution applies in some cases 
where such information is not available in the literature. 
We shall make use of the fact that a matrix of the form 


| 2 0 0] 
Qaub ante 


is positive definite if and only if a>0, b> 0, and bd>c?. Also, if wu; and v; 
are two vectors such that u;v; = u; v;, and if v; is not the null vector, there ~ 
exists a scalar A = u, v'/v! v; such that u; = A v;. If u; = u; and v; = v,;, where © 
u and v are scalars, then A and u must be functions of v. If uw; and v; are both © 
null vectors then, of course, u; = Av; holds for any A. 


2. Differential Equations Determining the Exceptional Tensors ¢;,; 


from 412) and: (1.3), 


Ap Dre ah as Pu 022) OD i. 
Os ae eam 1 J Sirti, |: te aA rs J 
z ar © i 2 rr Or 2 ors Ly hy par oe: 11,,) Coe 
fee 
0272 ( 
=2(srap ly tore U i) 4 | 
Differentiating (2.1), we obtain 
ii gaene Ono” Salty Woy a2 g* OA 
x Pik = Me Ge + a or + Me oye + Ge oparr + Se oree 
2. 2) 
038s 08s Ors’ 08a’ ( 
a’ Py ee ee 9 ; 
+ 4% ore + a orrory + Ge prone + Me SEE 
where 
ai, = (Ce er 
ai Cae , 
3 = = AR 
a;,=1, (¢ ae — as (c aA , 
Ay -1)j 
a;, = IL, (c Ae als Ci > [Il , (c aA a i ay k? 


Ee = oe a (IL; ah ) 
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Qe, = IT, L; (o> he I, U1, (c=) = 1721, él, 
ap UL (eS 1, Tl Ly IDs ot} 
a;,=11,11,¢. 


In order that there exist a function p satisfying (2.1), it is necessary and 
sufficient that the right member of (2.2) be a symmetric tensor. Suppose ¢;; 
is such that this condition is satisfied for every choice of the function (Ts ID). 
At a given point P, I and II have definite values, and we can choose »' such 
that all but one of the derivatives of Y occurring in (2.2) vanish at P. It then 
follows easily that the tensor which is multiplied by the non-zero derivative 
must be symmetric at P. It is thus clear each of the tensors a},, ..., a, must 
_be symmetric. We must thus determine symmetric tensors c;; which satisfy 
(1.4) and (1.9), which are such that Zh matrix ||c;;|| is eect definite, and 

such that each of the tensors a},, ..., a?, is symmetric. 

Suppose that I and II are not both constant, and exclude isolated points 
where I ; = II; = 0. The condition that a?, be symmetric is 


I, Linger LL cae 
and II ; = 0 implies that IT ; c, = 0. Thus there exists a scalar A such that 
II ; Ga Ail. (2.3) 


Either II ;= 0 or I; is an eigenvector of c} and A is the corresponding eigen- 
value. baie the bie avectors of c; and its inverse, (c~ 1)t, coincide either II ; = 0 
or II ; is an eigenvector of (c-¥)i. A similar analysis applied to a’, shows that 
either I ; = 0, or I; is an eigenvector of (c~1)?, hence of c}. Thus 


lg A l,, (2.4) 


where A is the eigenvalue corresponding to I ; unless I; = 0, in which case A 
is arbitrary. Setting the anti-synimetric part of a’, equal to zero, then using 
the above results, we obtain 


mets ilies nt Ea LP je (2.5) 


We assert that we can take A = A. This is clearly true if I; = 0 or I1,= 0, 
for then A or A is arbitrary. Suppose neither I,,; = 0 nor II; = 0 and that 
A + A. Then (2. 5) implies that I ; and IIL; are iatallal! But I ; Gnd II ;, being 
eigenvectors of ci corresponding re cae eigenvalues, are merpendiculae, 0) 
we are led toa PU aeatcnot Thus A = A. We next assert that I and II must 


and if all the eigenvalues of c} are distinct, then I; and II ;, being eigenvectors 
of ci corresponding to the same eigenvalue, are parallel. Thus I and II are 
functionally dependent. Let c,, cy, cs be the eigenvalues of cj. Then, from (1.7), | 


ill if ; 
ieee ae ays ey (2.6) 


C1 ls C2 C3 C3 Cy 
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be functionally dependent. If I or II is constant, the assertion follows. If not, 
| 


and, from (1.4), 
Cpeyty = T. (2.7) @ 
( 


In a region where c! has equal eigenvalues, it follows from these equations | 
that I and II are functionally dependent, so the assertion is true. Thus there © 
exists a nonconstant function B such that I = I(B), II = II(B). From (2.6) and i 
(2.7), we then have c, = ¢,(B), cy = ¢,(B), cs = ¢3(B). We assume the notation ~ 


chosen so that A = A = ¢,. Then (2.3) and (2.4) yield 
7 
Bec, Be. (2.8) 


Now 
a, = = B, Ge, LB O) x 
=—Il'B, On =(U' 2, Bo), 
Sec B cd = {Ll a) 0B, By ee a, 
where primes denote differentiation with respect to B. Thus, in order that a?, be . 


symmetric, we must have 


Il'(B , dh eg oh ch 3) =. 
Thus either II is constant, or there exists a scalar D such that 
Oe aD Be (2.9) 
A similar analysis applied to a?, yields 
V'[B, (coe eee (cr \hi] = 0, 
so either I is constant, or there exists a scalar E such that 


(cc) ,=EB,. (2.10) 


In order that aj, be symmetric, (2.9) and (2.10) must hold when I or II is 
constant, so these conditions hold generally. For the anti-symmetric part of 
a, to vanish, c} ; must be the gradient of a scalar. This condition, together 


with (2.8), implies that D = D(B). Similarly, the symmetry of a}, requires 
that E = E(B). 
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B. Now ci and (c~1)} can be written in the form 

, C= 6, a a, + 6, 01d, + cg! CA (2511) 
; (c~); ae = a’ a, + = bi b + = c C, (2,12) 


where a’, b', and c‘ are mutually orthogonal unit eigenvectors of ci. It follows 
from (2.8) that we may take 


Bis 


: a; ee ey Pe, 13 
VBE B, a) 

Also, since a’, b', and c’ are mutually orthogonal unit vectors, 
Gj =@a+¥b+cc. (2.14) 

From equations (2.11) to (2.14), 
Bute a 

Cj; = (Cy — C3) Be B - + (Cy — €3) b; b; + C3 8; (2.15) 
-1jij _ (1 i De tat A Gt 4 bi pi 1 oii 2.16 
(c ) a os | BEB, (= =) Sis ex g . ( . ) 


Substituting the expression for c;; given by (2.15) into (2.9) and using the fact 
that B ; and 6; are orthogonal, we obtain 


BiBi 
BEB, 


(ey — 0) (BoB *) + (eo — 9) (OB), = (D— 6) BY, 


while (2. 10) and (2.16) yield 
ee er ta cc) Os B-()e 


In order that these equations hold simultaneously, we must have 


(grg), 7 58) B ifcy+c3, (b'');=G(B) B’ iftg+cs, 
AB y |i 


where F and G are rational functions of ¢,, cy, ¢3, D, E, and ¢}. 
To recapitulate, we have shown that, if I and II are not both constant, the 


following conditions must be satisfied : 
¢,=¢,(B), ¢.=¢(B), ¢3=¢,(B), (2.17) 
6 Cote 1; (2.18) 
VS 1, (2.19) 
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(bb) , =G(B) Bi ifc, +6,, 


Also, |\¢;;| must be positive definite. If these conditions are satisfied, the 
tensors Ci and (c c~1)"1 given by (2.15) and (2.16) are such that each of the 
tensors a@},, ..., a), is symmetric. Furthermore, (1.4) is satisfied. 

If I and II are both constant, we can not be sure that these equations must 
hold. If they do, it follows from the analysis made in the following sections that _ 
c;; is a constant tensor. It seems very likely that a deformation for which I, II, . 
and |c}| are constant must necessarily be homogeneous, so that c;; is a constant 
tensor. Certainly one can always determine a homogeneous deformation for 
which I, II, and |ci| take on prescribed values, provided only that these 
values be consistent with the requirement that | c;;|| be positive definite. It 
would be of some interest to prove a general theorem from which the result © 
conjectured here follows as a special case, namely that the invariants I, II, 
and |c| determine the deformation to within a rigid motion of the material | 
as a whole. It is well known that, if these invariants are equal for two tensors 
¢ = |¢;;| and ¢ = |/¢;;|, then there exists an orthogonal matrix a such that 
c=a"'ca/[8]. The difficulty arises in showing that, if ¢ and € satisfy (1.9), 
than a@ is constant. 

Since solutions corresponding to deformations for which c;; is constant are 
discussed in [1] and [2], we shall not discuss them here. 


3. Remarks Concerning Solutions of Equations (2.19) to (2.22) 
From (2.13) and (2.21), 
aa +a,a'=FVB* B, a , 


so a’ and a’, a‘ are linearly dependent. But, since a’ is a unit vector, a’ and 
a’ ,a’ are perpendicular, hence linearly independent unless 


@ rae SO (3.1) 
Thus 


=F VB BS: (3.2) 
Equations (2.13) and (3.1) yield 


, 


(BY By 21 BeBe BeBe) es 
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“from which it follows that B* B , is a function of B,VB* Bh = H(B). Then 
. 2) becomes 


a,=F A. (3.3) 
Balso, 
4 Bi Bij H’ 

an = (3 ) = a BB; (Pa eae (3.4) 


_ Taking the divergence of (3.1), then using (3.3) and (3.4), we obtain 
f Poe) = es, AY '+a@ia ¥ a= Tt) Bae ata: sae. (3.5) 


We assert that the surfaces B= const are surfaces of constant mean and 
Gaussian curvature. The constancy of the mean curvature follows from (3.3) 
and the fact that + a',; is the mean curvature [9]. The Gaussian curvature 
is given by!) é . ha 
< [(a’,)?— a a,,—a', a a, a‘), (3.6) 
. 
“which, by (3.1), (3.3), and (3.5) is a function of B, so the assertion is proved. 
_ Suppose the mean and Gaussian curvatures of a surface S are constant. 
Then the principal curvatures 1/9, and 1/0, of S are constant?). Introduce 
surface coordinates (uw, u?) such that the coordinate curves are lines of curva- 
ture. Then, on S, the metric tensor g and second fundamental tensor d are of 


the form 


T 

|| G19 | | ne 0 
g= ‘ las d= " Son ) 

tI bee | Q2 


and the Mainardi-Codazzi relations become 


0 1 1 ( 1 1 Olog G11 __ 0 


Ou> 01 2 \@1 2 Ou? 
0 : 1 a 1 ( i 1 Olog £22 boi 
Our Qs 2 \@2 01 Our 


=) That (3. 6) gives the Gaussian curvature of a surface with unit normal a‘ can be seen by 
taking the equation of the surface in the form z—/f(x, y) = 0, where (x,y,z) are rectangular 
Cartesian coordinates, and using the fact that 


1 
p q : ie 
= - eee a: Se soe 
Titre ; yite+e * Vis p+ ¢ 
where p = 0f/0x, q = Of/0y. The expression (3. 6) then reduces to the well known form 
rt—s® 07+ 0? 07 


—————— = ? s= ==> ‘= ——. 
( + pet gy?’ where 7 Ox®, Ox Oy Oy? 


2) The equations and theorems of differential geometry used in this discussion may be found 
jn [10], pp. 150, 217 — 221, 225, 228 — 230. 
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Since 1/0, and 1/9, are constant, either 1/9, = 1/02, in which case S is a plane 
or a sphere, or 24, = £11(#1), Bo = &o9(#"). In the latter case, we can Tepara~ 
meterize the coordinate curves so that g1, = g9, = 1. Gauss’s equation oe 
requires that 1/9, 0, = 0. Assuming the notation chosen so that 1/e, = 0, we 

have 


1 | ; 

1 0 <= OH ; 

p84 =o , d =e Qt , 

Ora 0.0 : 

where 0; is a constant. It then follows that, in a suitably chosen rectangular ; 
Cartesian coordinate system (x, y, z), S is given parametrically by : 
1 ub 3 


U ie 
%=0,coS—, Y= o,Sin—, z=u?, : 
Qi Qi : 


so S is a right circular cylinder. Thus the surfaces B = const are spheres, 

planes, or right circular cylinders. 
Now consider the curves x‘ = x‘(s) satisfying the differential equation 

dx* 

ds 


=a. (3.7) 


Taking the absolute derivative of (3.7) and using (3.1), we obtain 


D LE RS er ee Fee ees 
dss eds eas ee Pe us 


which is a sufficient condition that these curves be straight lines. Thus there 
exists a congruence of straight lines normal to the surfaces B = const. For this 
condition to be compatible with the results obtained above, the surfaces 
B= const must be concentric spheres, parallel planes, or concentric circular 
cylinders. 

Let us digress for a moment and consider a problem in fluid dynamics 
which is solved by this analysis. If there are no body forces acting, the steady 
irrotational motion of an ideal fluid with a barotropic equation of state is 
governed by the equations 


Bay BY = Pie 0 Old) Ol Oe (3.8) 


where B is the velocity potential, g the density, and # the pressure. It follows 
from these equations that 


fess + = v? = const (3.9) 
and 
Bee): tae angie 2 4@ 3, tt 
(ar ave ys (2 Te 7 oe? nee (3.10) 
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where v = B* B, is the speed. Consider the problem of determining all flows 
of this type for which the speed is constant over the equipotential surfaces, 
v = v(B). Then (3.10) reduces to 


BiBi\ of Lap Wes 
ee Then eee 


But, from (3.8) and (3.9), o and # are functions of v, hence of B, so 


On SR ae \ ates 
(1+ 08 5) = F(B) 
Thus the equipotential surfaces are concentric spheres, parallel planes, or con- 
centric circular cylinders. Since these flows are well known, we shall not discuss 
them further. 


Let us now consider (2.19), (2.20), and (2.22). From (2.22), 
0 +0, b'=G(B) B’. (3.11) 
Multiplying (3.11) by 6;, summing on 7, then using (2.19) and (2.20), we obtain 


bY, =0, (3.12) 
so that 


ioe obs © 
b B= (6, — 8) B= (Si - Sa) = By, (3.13) 


where P = i G dB. But (3.12) and (3.13) are the equations governing the 
steady motion of a fluid of unit density, subject to no body force, if 6; be inter- 
preted as the velocity vector, P the pressure. By (2.19) the speed is constant. 
Suppose there exists a family of surfaces y = const to which }; is normal. Then, 
by Prim’s generalization [11] of a result due to Hamet [12], y = y(a0 + 62), 
where a and 0 are constants and (r, 0, z) are cylindrical coordinates. From 
(2.19), b; is a unit vector, so 


Va 
= 
Vvjnv* 
Thus ; 
ay Y 
= —_ = —— — 4 ———: Sete a 3h, 14 
= 0, b= La” b, aE (3.14) 
Equation (2.20) then yields 
0B OB 
tepy ab ety 
so that 
B= Bir,az—br* 6). (3.15) 


Since this result does not hold unless there exists a family of surfaces to which 
b; is normal, we shall not use it when a more general analysis can be made. 


4. The Case where the Eigenvalues of c;; are Distinct 


When the eigenvalues of ¢;; are all distinct, (2.21), (3.12), and (3.13) hold. — 
As was shown in section 3, the surfaces B = const are then concentric spheres, — 
parallel planes, or concentric circular cylinders. Suppose they are concentric 
spheres and choose a spherical coordinate system x1 = 7, 7 os Oa ag in 
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. 


which B = Bir). From (2.20), | 
| B,¥=B,Y=0,. 

so b’ = 0. From (2.19) 

g1 b, b, as He 
in, T 73 sin?0 
We thus have 
b,=0, by=rcosy,, b,=rsinOsiny,, (4.1) _ 
where y, is some function of the coordinates. Equations (3.12) and (4.11) yield ~ 
joae wl Osin 8 cos py Osin y, |: 

ie para | 00 eee |=0. (4.29 


From (3.13) and (4.1) 
RN es Ee a er 
(se — ait) or = a) = 0, 


Obs 5 0b, = [see Osin 6 anv] 2p 
a Op 06 tik 


0) 


(4.3) 


Now (4.2) and (4.3) require that 


cos | Lao, 


ows 
a 
sin py (cosd ae a) + COS, sino © =0. 


For these equations to be satisfied simultaneously, we must have 


Ow, 


pera OW, 
00 Ca. 


Op 


= —cos@. 


Since these equations admit no solution, this case cannot occur. 

We next consider the case where the surfaces B = const are parallel planes 
and introduce a rectangular Cartesian coordinate system x! = x, x2 = V), 08 ane 
in which B = B(x). By an analysis similar to that given above, we find that 


b,=0, b,=cosp,, b,=siny,, (4.4) 


> 
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where y, = (x). Using (2.15) to (2.18) and (4.4) to calculate ¢;, and (c~})‘/, 
_ we find that 


/ i | 

wy! lel = | 0 H J : (4.5) 

4 hee ee 

; HK—J? 0 0 

| che A K Ss 
fer L=) emo re |, 

| 0 aie ga 

: | HK—J* HK —f? 

: where 

- H = (Cy — Cs) cos* yy + C3, 

J = (Cz — €3) COS SiN Yo , 


K = (cz — cz) sin? y, + 3, 
are. functions of x. Making these substitutions in (1.9), we obtain, after some 
simplification, 
ee TT ae Ke fp a BK, 
where primes denote differentiation with respect to x. Thus 
H=ae4+a,, J=a%+@, K=a,x+a, 
H K — J? = (a, 4g + aq a5 — 2 dg a4) X + Ag Ae — GZ, 

where the a’s are constants such that 

az = a, a;. (4.6) 


We exclude the case where c;; is a constant tensor, that is where a, = a, =a, = 0. 
We now introduce a rectangular Cartesian coordinate system (x, y, z) such that 


x=x+1, y=yoosu+Zzsnu, z=—ysinu+Zcosp, 
where 


ai — a, A 
A, 4g + 4,45 — 243 a4 


, tanp==Z ifa,+0, w= > ifa=0. 


This transformation is not defined if 


Ay Ag + 4g 45 = 24344. (4.7) 


= io 
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4 
If (4.7) holds, then by (4.6) and (4.7) ; 
az az + 2 a3 (dy ag — a4) + a az =0. (4.8) 
- 
Since ||c;;|| must be positive definite, HK — J? > 0, so that, in this case, H 
Ay a, > aj. Thus each of the three terms in (4.8) is nonnegative and a, a + 0. — 
Thus a, = 4; = a;= 0, a possibility which we have excluded, so the trans- 7 
formation is defined. Transforming c,; and (c~1)"? to this new coordinate system, — : 
and dropping the bars on the variables, we obtain . | 
| 
by by x = e é 
Lot =i Tos. beaahe hy ae (4.9) 9 
0 bs b, . 
b, by x 0 0 

.; 1 —bs 
cat yy Seis ©2 oe SEA (4.10) 

0 —b, by % f- bs 

b, b, ¥ b, by * 


where the ’s are constants such that b, by = 63. In order that || ¢;;|| be positive — 
definite, it is necessary and sufficient that 


ee Out Oe, 
From (1. 7), (4.9), and (4. 10), 


isd b,x + by + by 2s 1 
T= 0), byx+ b, by x , Be ae Lk ak eh 


Equations (2. 1), (4. 9), and (4. 10) then yield 0f/0y = 0p/0z = 0 


? 


Op OZ dby byx oer di ey dil 
i deal DW eee Pipi tan ae STS as battin 
Roe ee ol ee OAs dl oxy di 1 
OLL. dx. 8b, b, 4 OI OIL dx * O11? dx b, b,% 


d [0s ox 1 
pes pe ees 
or 1% — ary mane: 


‘so 


az os 1 
Pie Yea See re 
Py Pot Sogn lathe varranopnee et: 
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where #, is an arbitrary constant. Using (1.3) to calculate the stresses, we 
obtain © 


tex =—Po> fie = 0; 


Ses 0x Ox 1 ox 
ey Por 4 a Oy + sr) Es b,x b, x) eee oie 
= oie f oz (4.11) 
tye = —2 0s (Fy sae ee aa 


0x 0x i b Ox 
Les =—pot 2 (5p e+ 52) Eee — by) +2 RE rae | 


Equations (4.9) and (4.11) constitute a general solution to the equations of 
this theory. In order to gain a better understanding of this solution, we shall 
determine a deformation compatible with (4.9). It is convenient to take the 
coordinates X* to be cylindrical, X1= R, X?= O, X= Z. Then, by (1.5) 
and (4.9), 


bees Maral lia) circa 


aR ee ee ZZ _ 


dx © dy dy ' Ox” dy 
Be ee 

(ae) +B (G5) + (Gp) het oe, 

Pe oy a ay oe eee, 
Fase Cele Car): 


These equations are satisfied if 


R=2]/ a O=b,V% 2, Z=+Vb,y+Vyz, 


b, b, ‘ 
from which 
Rr? ee pees 2 ee 4.12 
BAD ie wth emda rie aa: LM 


The orientation of the cylindrical system with respect to the Cartesian system 
can be assigned arbitrarily. Changing this orientation is mathematically equi- 
valent to subjecting the material to a rigid motion before deforming it. The 
deformation (4.12) can be described as follows: We start with a portion of a 
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right circular cylinder bounded by the surfaces 
Rim hy, k= Ky, 0=+090,, Z=+2Z,q. 


This solid is first deformed into a rectangular parallelepiped bounded by the 
planes 


R? Ri 0 
=o = dD, ri % tg =O", alias ag AS Se z= +Z,: . 


x-axis which deforms it into the rectangular parallelepiped 


7 


© =X, = OF Na Gea ais Km, | : 

} 
igi: S Vo é jeotigivr 5 (4.13) © 
Y= on Yor ser) Ac ence | 


so b, is the extension ratio, measured from the position occupied by the body — 
after the first deformation. Finally, it is subjected to a simple shearing defor- 
mation which deforms it into the solid bounded by the planes 


K=%, X=X%., YHHYo, 


and 
BS Mao = CSO 
r= t3,-/%y if b..> @, r= tat |» if 6, =<. 07 
4 4 


Thus Vb,/b4 is a measure of the amount of shear. 

We must now consider the case where the surfaces B = const are concentric 
circular cylinders. It is convenient to introduce a cylindrical coordinate system 
xt=y7, x? = 6, x= 2z in which B= Bir). By an analysis similar to that 
employed in discussing the case where the surfaces B = const are concentric 
spheres, we obtain 


b,=0, b=rcosy,, b,=siny,, (4. 14) 
where ys is a solution of the equations 


ee 0 0 
— sin ps ae + 7 COS Ps = =0, cosps fi + rsiny, ov =O: 


As is easily verified, these equations are satisfied if and only if ws = y,(r). 
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Using (2.15) to (2.18) and (4.14) to calculate c;; and (c~1)"7, we obtain 


i | ¢; | as 


I (e-*)""|| = 


where 


M = (c. — €g) 7 Sin ws COSY, , 


N = (€y — Cs) sin? ws ++ Ce , 


72 
LN—M 0 0 
0 L M 
0 M N 
LN—m: 
v2 0 0 
N —_M 
LN? sCW ow 
a SM Pe 
LN—M? LN— M2 


L =[(c, — €g) cos? + ca] 7? , 
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are functions of 7. Making these substitutions in (1.9), wo obtain, after some | 


simplification, 


Wis 
" 
LS = ee 


MM’ N’ 


Ui i] 
Mice Seay NY ceed, 


Y Y 


where primes denote differentiation with respect to 7. Thus 


Lil, "as 


LN — M? = (a, dg + Gq a, — 2 dz M4) 7? + Ay Ag — Gj, 


where the a’s are constants such that 


For I ¢;;|| to be positive definite, it is necessary and sufficient that 


af? +4579, 


Suppose first that 


Then from (4.16), 


ZAMP V/31 


az = 4, a,. 


(a, 4g + de ds — 243 a,) 7? + Aga, — az >0. 


Gg Gs\. 


(M’)2 = [,'’M' ‘ 


N =\2, 77. de 


(4.15) 


(4.16) 


(4.17) 


(4.18) 
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Squaring both sides of (4.17), then using (4.15), we obtain 


a? az + 2 (a, a4, — at) a2 + a az=O0. 


sere cane eee 


By (4.18), each term in this expression is non negative and 4, a, +0. Thus 
a, = 43 = 4; = 0. Furthermore, by (4.16), a, > 0, a, > 0. Thus there exist real 


constants 0,, 6,, bs such that 


Eee aot a a 
b, = V ay ag — AF b, = |/ ——_* = ba = —4; 
ii Vay a 4, 9% ip, ee ae 


It then follows that 


oe 0 0 
lesl=|) o 24 H Ps Bb, | (4.19) 
0 BRB, 2 Be 
ae 0 0 
| (c-*)*"|| = || 0 b5 —b3b, ||, (4.20) 
0 sit 1 tae bs 
le ee ee r= Semen re 


Equation (2.1) then yields 04/00 = 0p/0z = 0, 


op _ 70k a of a (ae dl , 0o8F aI) b3 
Or Ol dr 7? Olt Aire al oul > on 


BDO zt Cake amp dl oy dil) v2 
oll dr BB Ol DILL era Olle = b2 


y2 


Damen OF kee oy 1 rt 14 d9bgd3 

ar aes A 03) a ia Cae & aera #1 ‘) 

isan DAU? haa a 

dy [2 et ant b3 Z|, 
so 

O27 bt Os 93 
= | ae 
DS Poth Are aa aera b} 2 ay 


where #, is a constant. Using (1.3), (4.19), (4.20), and (4.21), we can compute 
the stresses. The physical components of the stress tensor are found to be 


| 
: 
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given by 


ee — boat f= 0 72 = 0; 


0-292 (ye) 2 ELEM 
=a (Ep), (4.22) 
Bm oe (apg et) 2 one (ap — 84) +2— de, 


If we choose the coordinates X% to be rectangular Cartesian, X1= X, X2= Y, 
X* = Z, we may write a deformation corresponding to (4.19) in the form 


Zz 


ra V2b,X; 6=3,Y, Sa 
i Diet| 


This deformation can be decomposed into two parts: we start with a rectan- 
gular parallelopiped bounded by the planes X = X,, X = X,, Y=+ Yj, 
Z=-+ Z,. This solid is then deformed into the solid bounded by the circular 
cylinders e 
| pep V2 BER Pg V 2b, 


and the planes 


Zi 
6=+O=145%, este Se 


A solution corresponding to this part of the deformation is discussed in some 
detail in [4], where it is shown that, if the bending is to be accomplished by 
terminal couples, the constants 6,, b,, and fy should be chosen so that 


b= tb. 7%, Pyo= (1). (4.23) 


In the present case, the bending deformation described above is followed by a 
deformation which carries the body into the solid bounded by the cylindrical 
surfaces r=7, and r=7,, the planes 6 =-+ 6, and the helical surfaces 
z+ bs 0 = + 2. By (4.22), the surfaces y = 7, and ry = 7, are subject to normal 
stresses only, the planes 0 = + 6, to a shearing stress in the direction of the 
z-axis and a normal stress, the helical surfaces to a normal stress and a shearing 
stress in the direction of the helices obtained by intersecting these surfaces with 
the cylinders y = const. One can always choose fp so that the surface 7 = 7 
is free of stress. Suppose this is done and that 0, and }, are chosen in accordance 
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with (4.23), so that I(7,) = I(7,). Consider a neo-Hookean material, that is\4 
material for which X' = « (I — 3), where « is a constant. In this case the surfac 
v =r, is free of stress and the resultant normal force acting on the plan 


06=+ 4, 


N = 22% [ 60 dr = 2 zr (Z — Py) se 
' 
vanishes. For a general form of the strain energy, these conditions are not 
satisfied. In general, these will be normal forces acting on these surfaces. If 
these are not supplied, the material will tend to expand or contract in these 
directions when the latter part of the deformation is effected. : 

Suppose now that (4.17) does not hold. If a, (a, ag + a, a; — 2 a3 a) < 0, 
then from (4.15) and (4.16), 


4 


(a, dg + a2 a; — 2 Az ay) (ft eS) (a, ag + Ag a, — 2 ag ay) 


2ids Ane @2 he 
eae 2 2 2 43 M4 345 
= Ay — Ap Mg — + Sie one 
ay ay 

a aza 

S aj — a, a, — at + 2 |a, a,| //— — == 

4 2 U6 4 2%, a ‘a 

1 1 


= 05 ~ 054 —|Jo4] — la] |/2 | 


9 
S aj — aga, 


which contradicts (4.16). Thus a, (a, a+ a, a;—2a3a,) 20. Similarly, 
As (ay ag + Ay 4; —- 2 az ay) = 0. Using these conditions, it is easy to show that 


there exist real constant 0,, ..., bg such that 
by = 4, Ag+ 4,4,—2a3a,, b,b3b,= —az, 
by by = a} — ag Ag, Db, be by by = 4, , 
b, bs = a, 6, 3 =a,;, 
bf 07 — by by 62 = ay), bi bg — b, b, 62 = a, . 


In terms of these, 


Wh | 
by V2 — by) : : | 
0 b, [by 62 + b2 (7? — b,)] — by [by by be ++ by bs (x? — b,)] 

0 —b, [b, by bg + b, bs (7? — b,)] by [b, bf +- b2 (7? — be) | 
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The solution corresponding to these strains is discussed in [6], where it is shown 
that deformations corresponding to (4. 24) are possible not only in all isotropic 
incompressible materials, but in all incompressible materials possessing trans- 
verse isotropy.*) 


5. The Case where c, = C; 


If cy = cg and ¢, = cs, then by (2.18), cy = cy = cz = 1, so the material is 
unstrained. If c, = ¢3, but ¢c,-+ cy, (2.21) holds, so the results obtained in 
section 3 apply. Thus the surfaces B = const are parallel planes, concentric 
circular cylinders or concentric spheres. Furthermore, from (2.15) and (2.16), 


Bor (5.1) 


1 1 Bt Bi 1 


BB; -1)ij 
a (Cy — ¢3) 7B ce ie 6s) BEB, " 


B 
where © 
OE oe be (5.2) 


Suppose first that the surfaces B = const are parallel planes and introduce a 
rectangular Cartesian coordinate system «x! = x, x*=y,x%=z in which 
B= B(x). Then from (2.17), (5.1), and (5.2), 


mo 
ej] = | One O ' 
| O08 Cy ol 


where Cs = C,(x). Since this is a special case of (4.5), we obtain no new solutions 
here. Similarly, no new solutions are obtained in the case where the surfaces 
B=const are concentric circular cylinders. Suppose they are concentric 
spheres and introduce a spherical coordinate system x4 = 7, x? = G, 4° =@ 
in which B = B(r). Then (2.17), (5.1), and (5.2) yield 


ae, Gr sii | @ 0 0 
e® 1 | 
= INERT | 1 
Le; || = | O--co#? 0 ley | = | 0 ee 0 P 
: i 1 
| Ones 0337 'sin* 0 | 9 0 art sino 


where cz = ¢3(7). Making these substitutions in (1.9), we obtain 


€3[(¢3 77)? = 477, 
4) Note added in Proof: The bending solutions discussed in [13] are included among those 
tiven by (4.22) and (4.24). 


. 
tq 
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from which cz = 7-2 (r? + a), where a is a constant. Taking the coordinates” 
X* to be spherical, X1= R, X?= O, X?=@, we may then write a correspond-_ 
ing deformation in the form : 


y=(-R?— a)", 0=+0, o=O. 


A solution corresponding to one of these deformations is discussed in [5]. 


6. The Case where c, = Cc; 


Since there is no essential difference between the cases c, = cz and c, = C5, 
we shall consider only the former, assuming that c, + c,. Then, from (2.15), 
(2.16), and (2.18), 


1 1 


as 186 uf 
C;; a (Cy == Ca) b; b; = Cy Si; ? (ery? — (— a =) oY o cs 


—gi, (6.1) 
5 


where 
Coste aes Ae (6.2) 


We restrict our attention to the case where there exists a family of surfaces 
to which 0; is normal so that (3.14) and (3.15) hold. Let us first consider the 
case where b = 0, so that B = B(r, 2). From (2.17), (3.14), (6.1), and (6.2), we 
obtain 


e* 0 0 
t=O ee 071), (6.3) 
0 0 e 
ee 0 0 
[ery] = es | 
Ome 0. eee 


where A = A(r, z) = logc;. From the conditions of compatibility (1.9), we obtain 


02a Oa By on 
Oe ot 3 
~ Or Oz (5 3 >) Oz’ (6.4) 
024 hy a 0 x 
Ors et (6.5) 


oar =3(e) -2. 4 (Sy), (6.6) 


(Bet Be -eP on 
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In order that these equations hold simultaneously, 2 must be constant. If 
04/0z = 0, this follows immediately from (6.5) and (6.7). Suppose 0A/dz + 0. 
Adding (6.6) and (6.7) and using (6.5), we obtain 


an fete oa 6.8) 


Or Oy, o@ oy ee Or 
Differentiating (6.8) with respect to z, then using (6.4), (6.7), and (6.8), we 
obtain 


dil OA rene ee (52) 0, 


y Or v2 Or 


(6.9) 


from which 0?4/(0r 0z) = 0. But from (6.4), we must then have 04/dr = —1/r, 
which does not satisfy (6.9). Thus 4 must be constant. It then follows from 
(6. 3) that c;,; is a constant tensor. 

_ Next assume that a = 0. We may then introduce a rectangular Cartesian 
coordinate system (x, y, z) in which 


b,=0, b,=0, b,=1, B=B(x,y). 


From (2.17), (6.1), and (6.2) we then obtain 


; | 
= 0 Of Pi 0x 60 
1 avs 
lei] =|} 0 a Ihe Cape cal OS ars a 
0 0 P Om Ot ae 


where P = P(x, y) = 1/cs. Setting R,,,, = 0 in (1.9), we obtain 
02P OP \2 OP \2 
Bea) ae): 
while R,,,, = 0 requires that 


ae er ee 
“Oy? ( Ox oy 


24 


Using these conditions and the equation R,,,, = 0, one can easily show that P 


must be constant. Thus c;; is a constant tensor. . 
Suppose now that a +0, 6 +0, and introduce a curvilinear coordinate 


system (&, 77, ¢) such that 


ee heir ate BEC 
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In this coordinate system, B = B(n, ¢) , 


Hees x : 
lal=|} 0 ayiee 2b6 eer | 
ic ROY M2 Gea pecs io Nes oe : 
| a i cP 0 ; 0 
i= | 0 temen fives Git) aac Sse 
0 200-827. 1 
b= aaa Soe b=, =0, Bap VEEP’ pivao. 


From (2.17), (6.1) and (6.2), we then obtain 


QO" 8 0 ‘ eae 
; ee ee 
Feab= 9 7 Cre CBnt |, ei) 2G ae 
gnt | 
0 Og “Oe | °c ome 
where Q = Q(y, ¢) = cg. Setting Rz,,<, = 0, we obtain, after some simplification 
AeOeig © OD eet Sa Oe 
Tate Bt 7G Oy Be (a& — bn)? | 
zs ~ lig, POs eae e. Zs 
Seas on) [o On Of  Q® (a? + b2 C2) 
b2 f3 0Q-2 


a®+b?f2 On (6.10) 


< btacss p oD a OG} 
Ineis AW gp oh 
4a2b2027Q 1 OQ (a®+b202) 0 C2 
(a? + bt 04)? — Q ac OC tats oeey | 


Now a& — by, n, and ¢ are functionally independent and Q = Q(n, ¢). Thus, 
in order that (6.10) hold identically in these quantities, the coefficient of each 
power of a € — b y must vanish, from which 


90° 9. epee? 6%) fae . 4.a2 b? 62 Q2 


0c ies Q? (a2 + b? C2) (a? + b2 £2)2 
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It then follows that R:;,: +0 unless Q = 1, in which case c, = c,=c,;=1, 
so the material is unstrained. 
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Résumé 


Nous avons dérivé trois solutions des équations de l’élasticité non linéaire. 
Nous avons déterminé les déformations qui sont possibles dans un solide isotrope 
st incompressible quelconque. 
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Die Konstruktion der mittleren Linie bei Oberflichenprofilen’ 
und einige Zusammenhinge zwischen den zur Kennzeichnung_ 
von Oberflachen benutzten Grossen 


Von Max Gary, Braunschweig?) 


1. Ubersicht 


Das Ziel der Priifung technischer Oberflachen besteht darin, diesen eine 
oder mehrere Kennzahlen zuzuordnen, die geeignet sind, diese Oberflachen zu 
bewerten. Es liegt dasselbe Problem vor wie bei der Untersuchung statistischer 
Massen. Es handelt sich um eine Mannigfaltigkeit von co? Gréssen, denen eine 
oder mehrere Zahlen zugeordnet werden sollen. Dadurch erwachsen eine Reihe 
von Fragen: 

In welcher Weise kénnen diese Kennzahlen eindeutig gefunden werden ? 

Welche Oberflachen werden umgekehrt durch diese Kennzahlen beschrie- 
ben? Das heisst, haben die erhaltenen Kennzahlen iiberhaupt einen Sinn? 

Und schliesslich, wenn mehrere Kennzahlen ermittelt wurden: Wie kénnen 
die einen aus den anderen berechnet werden? 


2. Profilschnitte und Flachenkennzeichnung 


Legt man durch eine Flache einen ebenen Profilschnitt (siehe Figur 1), so 
ergibt sich durch die Form der entstehenden Kurve ein Kennzeichen fiir die 
Flache. Die Flache lege dabei in einem rechtwinkligen Koordinatensystem, 
und die Ebene des Profilschnitts wie auch diejenigen weiterer Profilschnitte 
sollen senkrecht zu einer im folgenden als Basisebene bezeichneten Koordina- 
tenebene sein, namlich jener, iiber der die Flache ausgebreitet ist. 

Nun ist die Mannigfaltigkeit der Geraden einer Ebene oo?, und auf einer 
jeden dieser Geraden gibt es co! mégliche Lagen einer Bezugsstrecke S,, so 
dass also insgesamt oo verschiedene Profilschnitte der Lange S, herangezogen 
werden kénnten, um eine Aussage iiber die Flache zu treffen. Diese Profil- 
schnitte schneiden sich zum Teil oder fallen gar teilweise zusammen. Es ist 
offensichtlich, dass zuniachst nur solche Profilschnitte, die héchstens einen 


Punkt gemeinsam haben, zur Berechnung einer oder mehrerer Flachenkennzah- 
len benutzt werden diirfen. 


1) Mitteilung aus der Physikalisch-Technischen Bundesanstalt. 
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Zweckmassig scheint es zu sein, aus einer Auswahl von verschieden ge- 
legten Profilausschnitten, die sich selbst nicht schneiden, durch harmonische 
‘Analyse die Fourier-Koeffizienten zu bestimmen und, vorausgesetzt, dass deren 
Streuung ein noch festzusetzendes Mass nicht tiberschreitet, durch Mittelwertbil- 
dung oder durch die sonstigen in der Statistik iiblichen Verfahren neue zuge- 
ordnete Fourier-Koeffizienten zu berechnen, die als Kennzahlen fiir die ge- 
samte Oberflache angesehen werden kénnen. 


Figur 1 


Ebener Profilschnitt durch eine Flache. 


Sofern sich zeigt, dass bei den einzelnen technischen Oberflachenbearbei- 
tungsverfahren eine erkennbare Verteilung der Fourier-Koeffizienten auftritt, 
die in den Einzelfallen lediglich durch einen konstanten Faktor voneinander 
verschieden sind, ist es sinnvoll, diesen Fourier-Koeffizienten selbst wieder 
neue Zahlen zuzuordnen: 

Wenn die Fourier-Koeffizienten mit 


iD eared ee, 3, 2. 5 9) 


bezeichnet werden, entsprechend einer Reihenentwicklung — die Darstellbarkeit 
mit endlich vielen Reihengliedern vorausgesetzt — 


ie 256 : 22 
e s cae ] 1 
y= : (a, cosr 3, x + b, sinr S, x) + dy, (1) 
so lasst sich der Wert 


Agus = 


haus = VES E+ +4 a 


finden. 
Es ist dies nicht die einzige Méglichkeit, die Gréssen 


Ged, a= Le 2a 
eimey Zahl zuzuordnen. 
Bildet man zum Beispiel 
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in der Form 
: (3), 
; 


nN 


Poe lhl dedeay (4). 


: 
und lasst zu, dass die Harmonischen beliebige Phasenverschiebung haben, so” 
erhdlt man einen Wert R,,, der den maximal méglichen Abstand eines oder 
mehrerer Punkte der untersuchten Flache von der Basisebene darstellt. Be-_ 
zieht man hingegen die einzelnen Profilschnitte, aus denen die Fourier-Koeffi- 
zienten berechnet wurden, auf die «mittlere Linie», so ist R,, als «mittlere 
Flachenrauhtiefe» zu bezeichnen. Hierauf wird im folgenden Abschnitt noch 
naher eingegangen. 

Eine weitere Méglichkeit, den Grdéssen 


ey GreOnng (f= eel coe 


eine Zahl zuzuordnen, ist durch die Formel 


n 


Dy (a, cosr 7 x + b,sinr fab x) + a 
flaunt: Sp Sp 


” 


1 
have = gsi / 
0 


gegeben. 

Nachdem man sich durch Priifung der Streuung der Fourier-Koeffizienten 
Klarheit dariiber verschafft hat, in welchen Fallen es sinnvoll ist, eben diese 
Koeffizienten oder gar die aus ihnen berechneten Gréssen als Reprdsentanten 
einer ganzen Flache zu betrachten, kann man, wo dies also angangig ist, auch 
gleich aus den einzelnen Profilkurven die Gréssen hpys, h4yp oder R ermitteln 
und deren Mittelwert oder nach einer in der Statistik tiblichen anderen Methode 
einen zugeordneten Wert bestimmen. 

Hier tritt ein neues Problem auf. 

Wie eingangs erwahnt, ist die zu untersuchende Flache auf die Basisebene 
bezogen. Diese Ebene ist in Wirklichkeit bei einer vorgegebenen Oberfliche 
nicht realisiert. Sie muss erst eingefithrt werden, und zwar bezogen auf die 
Flache selbst. 

Nun muss hier eingefiigt werden, dass die Oberflache in der Praxis durch die 
aufgezeichneten Profilkurven, die selbst keine Bezugslinie haben, erst gegeben 


dx (5) 
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wird. Es muss daher fiir jede Profilkurve eine Bezugslinie, «mittlere Linie» 
genannt, bestimmt werden. Dass diese mittleren Linien nicht exakt nur einer 
einzigen Ebene — der Basisebene — angehéren, ist dann unvermeidlich. 


3. Bestimmung der mittleren Linie 


Im DIN-Blatt 4762 (Februar 1952) wird als Definition der « mittleren Linie 
L,,» angegeben: « Die innerhalb der Bezugsstrecke so gelegte Linie, dass die 
Summe der werkstofferfiillten Flachen tiber ihr und der werkstofffreien Flachen 
unter ihr gleich und ein Minimum sind». Ob diese Definition hinreichend ist, 
um eindeutig die « mittlere Linie » zu finden und nach welchem Verfahren dann 
diese Linie auch wirklich gefunden werden kann, ist bisher nicht bekannt. 

Bei dieser Definition erhebt sich namlich sofort die Frage: Soll iiber einer 
bestimmten Bezugsstrecke durch zwei Senkrechte in den Endpunkten der Be- 
zugsstrecke aus der Profilkurve ein bestimmtes Stiick ausgeschnitten werden 
und jetzt die «mittlere Linie» so bestimmt werden, dass die erwahnte Forde- 
rung erfiillt ist ? Oder ist die Profilkurve lediglich durch zwei Punkte begrenzt 
und soll die «mittlere Linie» dann entsprechend der aufgestellten Forderung 
gesucht werden? Im letzteren Falle ware die Forderung: « Die innerhalb der 
Bezugsstrecke so gelegte Linie, dass...» usw. unverstandlich. im folgenden wird 
ein Naherungsverfahren beschrieben, nach dem sich diese Aufgabe lésen asst. 
Allerdings ist tiber die genaue Lange der Bezugsstrecke zundchst keinerlei 
Aussage gemacht. Vielmehr wird lediglich verlangt, dass auf der Profilkurve 
zwei Punkte das zu betrachtende Stiick derselben begrenzen. Die Lange der 
Bezugsstrecke ergibt sich erst hinterher. Zundchst sei die Profilkurve (siehe 
Figur 1) in bezug auf das willkiirliche rechtwinklige Koordinatensystem (4, Vv) 
gegeben. Dann lasst sich durch Planimetrieren die Gerade L,,, finden: 


cS 
1 7 
vii ‘ie a oe / V1 (%4) dx, . (6) 
6 


Eine weitere Planimetrierung der Flachenstiicke, die auf einer Seite der Linie 
Lm liegen und die von dieser Linie und der Profilkurve sowie gegebenenfalls 
noch vom Lot durch den Punkt P auf die x,-Achse begrenzt werden, liefert den 
Wert F, (siehe gestrichelte Flache in Figur 1). 

Dreht man jetzt das Koordinatensystem um den Winkel «, der im Gegen- 
sinn des Uhrzeigers positiv gewahlt sei, um den Punkt O, dann erhalt man das 
neue Koordinatensystem (9, yp). Die senkrechte Projektion des Punktes P 
auf die x,-Achse liefert den Punkt Ag. 

Durch die Strecke OA, ist die Bezugsstrecke S, gegeben. Wiederum erhalt 
man durch zwei Planimetrierungen den /, entsprechenden Wert /, und die 
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neue mittlere Linie L,,,. Die Fortsetzung des Verfahrens, wobei der Koordi- 
natenanfang stets beibehalten wird, liefert durch Interpolation aus den zuge- 
ordneten Werten 
ay — 0 | he | hes | 
PR ee at 


3 
den zum Winkel «* gehérenden Minimalwert F*. Diese Interpolation wird | 
zweckmassig graphisch durchgefiihrt. 
Der Wert F, muss in jedem Falle ermittelt vee um den Maximalwert ; 
F* wenigstens genahert zu bestimmen. - 
Die graphische Interpolation kann auch durch die numerische ersetzt : 
werden: 


ot) (a5 = 4) (7) 


(a — 
Fy (a5 = ey) ey {ay = a3) ] 


Whee (F; — Fa) 
2 TAL Pog aah 


Der Wert F* folgt dann zur Kontrolle aus zwei weiteren Planimetrierungen. — 
Zur naherungsweisen Bestimmung der mittleren Linie sind somit minde-— 
stens sechs Planimetrierungen erforderlich. Die Kontrolle bedingt zusatzlich — 
zwei weitere Planimetrierungen. 


An Stelle der Definition der mittleren Linie entsprechend dem Normblatt 
ware auch eine andere Festlegung denkbar: ; 

Mittlere Linie = Achse des kleinsten Tragheitsmoments der Flache, gebildet 
aus der Profilkurve und der Verbindungslinie der Begrenzungspunkte derselben. 

Der Vorteil: Die Definition ist eindeutig. Durch zwei Planimetrierungen mit 
einem Quadratplanimeter erhalt man eindeutig den Schwerpunkt der Flache. 
Durch drei weitere Planimetrierungen mit einem Potenzplanimeter zur Bestim- 
mung von f y® dx erhalt man die gesuchte Linie. Samtliche Planimetrierungen 
kénnen mit einem einzigen Instrument vorgenommen werden}). 


4. Der Zusammenhang von hyyr mit Apys in Gestalt 
einer Abschatzung 


Bezieht man die Darstellung der Profilkurve auf die mittlere Linie (Gerade), 
dann verschwindet das absolute Glied in Gleichung (1): a) = 0. Des weiteren 
sei verlangt, dass nur eine endliche Anzahl von Reihengliedern auftrete. 

Mit Hilfe der Schwarzschen Ungleichheit 


Eas) Ea) h9 


1) Vel. L.A. Ort, Systematische Entwicklung der Planimeter und I ntegrimeter aus der einfachsten 
‘Grundform, Messtechnik 13, 41—48 (1937). 
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ergibt sich aus Gleichung (5) 


s — aa a Sane tt 
Ngve Ss | | yz (az + | by (cos*r “2 x + sin2r Se ‘| | es 18} 


1 


t é as i — 
1 we — 
ave s = ie yeu (a2 + 0?) Vn dx 
0 | 


it 


have S <Vn x (a? + 67) (9) 
Mit Benutzung von Gleichung (3), in der ay = 0 ist, folgt 


= (10) 


nm = maximale Anzahl der fiir die Darstellung der Profilkurve bendétigten Ko- 
sinus- oder Sinusfunktionen. 

Eine untere Schranke fiir 44), lasst sich fiir spezielle Falle von Profilkurven 
angeben. Unter der Voraussetzung, dass sich die Profilkurve in der Form 


: 2% 3 2% 
y= 4,sinr—-~ > X + a, sinus % (11) 
darstellen lasst, gilt 


bore 5 [rte 


wie unmittelbar einzusehen ist. Folglich: 


a, Sin ae x ax, (12) 
: I 


A, "la;siny 22 x Oi 


1 
Sy 
0 


have = (|a,| —|a Ae (13) 


Wenn es gelingt, fiir eine nach bestimmten Herstellungsverfahren erzeugte 
Flache die obigen Voraussetzungen (zwei Sinuswellen) als erfiillt zu erkennen, 
so kann fiir den Quotienten Apys/Aaye eine obere Schranke angegeben werden. 

Ist 
eat on > 1, (14) 

m 


so folgt 


und nach Gleichung (3) 
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, mea | 
Daher: jae ne ae ES ere (15 


Beispiel: m = 10. hems = have i 1,24. 
In der Elektrotechnik heisst die Grésse 


h E 
Eid Ee Formfaktor» = pat 4 


hAVE Emitter 
Einen Uberblick iiber die Werte hrys und hayeg fiir einige Kurvenformen mit 
den Bezeichnungen der Elektrotechnik gibt das Archiv fiir Technisches sae 
V3620-2. | 

Summary ] 


Attention is directed first to the analogy between the problem of evaluating 
collectives and that of finding a numerical expression for an existing surface. 

For the value heys/haye, approximations are made using Fourier series of 
the surface profil based on the center line. h4yg = center line average (= /cza).— 
For hayg, an upper bound calculated from the Fourier coefficients is stated. 
With this, and as hpys may be determined in a definite form, a lower bound for | 
hrms/have is obtained. Provided that the surface profil is the superimposition — 
of two harmonics, a lower bound for h4yg and then an upper bound for hrys/h4yve 
could be found. In order to get a base line for which the areas embraced by the 
profile above and below the line are equal and a minimum, a graphical method is 
constructed. 


(Eingegangen: 25. November 1953.) 


Anwendungen des Quotienten-Differenzen-Algorithmus 


Von HeErInz RUTISHAUSER, Ziirich!) 


1. Einleitung 


Nach einer ersten Mitteilung?) tiber die theoretischen Grundlagen des 
QD-Algorithmus sollen nun einige Anwendungen behandelt werden. Vorerst 
wiederholen wir aber die Definition des QD-Algorithmus: 

Mit den Ausgangswerten s, =-s” wende man nacheinander fiir vy = 0, 7, 2,. 
und dann fiir o = 1,2,3,... unter Beachtung von ey) = = 0 die folgenden Formeln an: 


Ss 
(eee os (ye 1 ) 2 
dren gt HP) eed Lehto 


o 


dy Tasttut fiir angewandte Mathematik der ETH. 

*) H. Rurisuauser, Der Quotienten- -Differenzen-Algorithmus, ZAMP 5, 233 (1954). Diese Arbeit 
wird im folgenden immer mit I zitiert; sie enthalt auch das gemeinsame Literaturverzeichnis fiir 
beide Arbeiten. Eine weitere Arbeit tiber Anwendungen des QD-Algorithmus (Bestimmung der 
Eigenwerte und Eigenvektoren einer Matrix) wird in Bd. 6, 1955, dieser Zeitschrift erscheinen. 
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Die so erhaltenen Gréssen werden wie folgt in ein Schema eingeordnet (QD- 
Schema der s,): 


5 
0 
a 
s) dq) = ef) s0) 
1 0 
g® g? 
s®) a — @)|| 5) do ef |) 5 
(2 1 
. gq qs” qs 
3 Zh 2 2 1 1 1 
O dal @ a | 
if 


Es sei an dieser Stelle daran erinnert, dass man fiir die Anwendungen haufig 
nur die Werte g” und e®) benétigt, man kann dann die Kolonnen s¥), df), s’?, 
d® ... weglassen und erhalt so das QD-Schema in komprimierter Form [hierzu 
die Formeln (I, 4)]: 


0 
s© 
0 
n 
3) 0) 
1 (0 
gh , qs : 
(2) of 200) 
2 a (0| 
q ?) GP : q3 
<3) e) ef) 


2. Umwandlung einer Potenzreihe in einen Kettenbruch 


Ist eine im Unendlichen regulare Funktion 
= Ss 
er 1) 


gegeben, so kann man aus dem QD-Schema der Koeffizienten s, nach (I, 9) 
unmittelbar die fiir die Aufstellung des S-Kettenbruches notwendigen Gréssen 


ablesen ; es ist namlich 


| | | @] 2] 


(2) 


ZAMP V/32 
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Das gilt in vielen Fallen auch dann noch, wenn die Reihe fiir f(z) tiberall diver- 
giert, beispielsweise fiir 


co 
ee ae ee) MR ha 
1(2) = a eae? Cie ras oe (sper Reihe fiir e af . dz). 
z 
Das QD-Schema lautet naimlich in diesem Fall: 
Sy Gm ea re a 
| 1 
-1 
-1 —1 ‘ 
oe = Allgemein: 
2 —1 —2 
. —3 ; =e! —3 gq? =—( + ») 
ERP a Sao be foe ber 
24 Fi —2 


Somit nach (2): 
(2) 


ey 


< =i | Z| 3) 3 | 
oc Maal at 


a te A a 


Dieser Kettenbruch konvergiert fiir alle z, die nicht auf der negativ-reellen — 


Achse liegen (vgl. [8], § 93). 
Soll eine im Nullpunkt regulare Funktion 


f(a) = D's, 2 
0 


in einen Kettenbruch entwickelt werden, so betrachtet man die Hilfsfunktion — 


ie poems oe 
uw ba) Say ec 


deren Kettenbruch nach (1) durch 


H 
i 
: 
i 
i 
: 


(2) fl P| Pl a 
ch ANSLL | w 1 | u | 1 
gegeben ist. Setzt man hier fiir w wieder 1/z ein, so fiihrt eine Aquivalenztrans- 
formation, nadmlich Erweitern aller ungeraden Teilbriiche mit z, unmittelbar 
zum gewiinschten Resultat: 

Se =e S| z qe) | z e| Zz q | zg 2) | Zz qo 

OF 5,8 el TR Re ee ae 
7 1 laa 1 b= La 1 
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Stellt man beispielsweise fiir die Koeffizienten der Funktion 
eet gz a 
Cee Nal i 3r at 


das QD-Schema auf, so erhalt man 


gy”) = i : g?) = Disp Oi : Aves er 
y+ a (y+ 20— 2) (v+2oe-—-1)’ @ (» +26 —1) (v+ 20) 


und damit durch Einsetzen in (3) den bekannten Kettenbruch 


Zz 
é ne 


Dies Z7| 2/2 | 2/6 | 2/6 | 2/10| 2/10| 
j1 | 1 fed eee < Rade 


3. Summation schlecht konvergenter Reihen 


Fiir viele Zwecke der angewandten Mathematik sind auch geometrische 
Reihen nicht gentigend gut konvergent. Beispielsweise bendtigt man zirka 


20 Glieder der Reihe 3’ 2~”, um eine relative Genauigkeit von 10-® zu erreichen. 
: 0 


Fiir unendliche Reihen, die sich ungefahr wie eine geometrische Reihe verhalten 
(das heisst, der Quotient von zwei aufeinanderfolgenden Gliedern variiert lang- 
sam), erweist sich der QD-Algorithmus als praktisch brauchbares Summie- 


rungsverfahren!). 
Dessen Prinzip stammt von STIELTJES®), der es insbesondere zur Summa- 


tion divergenter Reihen verwendete: 


co 
Ist Ss, die zu summierende Reihe, so entwickle man die Funktion 
0 


co 


/(2) ae Pes 


in einen J- oder S-Kettenbruch und setze dann z = 1 ein. 
Als Beispiel betrachten wir eine Reihe, die im Verlaufe einer grésseren 


Berechnung auftrat: 
0,775 551 + 0,587 903 + 0,451 730 + 0,350836 + 0,274795 + 0,216681 + ---. 


1) Wenn jedoch die zu summierende Reihe nicht den Charakter einer geometrischen Reihe hat 
(zum Beispiel 1 — 1!-++ 2! — 3! ++ 4! — +- ---), soist wohl das Verfahren von D. SHANKS vorzuzie- 
hen. (Vgl. D. SHanxs, An Analogy Between Transients and M athematical Sequences..., Naval 


Ordnance Laboratory Memorandum 9994.) ; 
2) SrieLTjes, (Euvres, Bd. 2; siehe auch H. WALL [8], § 93. 
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Wir stellen das QD-Schema auf: 


(v) (v) (») 
Sy gy? e? I ee q3 
0,775551 
0,7580456 
0,587 903 0,010 3295 
0,768 3751 0,6155258 
0,451730 0,008 2747 0,0129556 
0,7766498 0,620 2067 0,4790043 
0,350 836 0,0066079 0,0100060 
0,783 2577 0,6236048 ~ 
0,274795 0,005 2610 
0,788 5187 
0,216681 
Hieraus nach (2), wenn man z = 1 einsetzt: 
<S Ween | 0,758 0456 | 0,0103295 | 
Po at ee ieee vers 
0 1 
DVAS) 129556 0,479 004 3 | 
a ate ; | Oe ws 3,52348, 


was mit allen angegebenen Stellen mit dem auf andere Weise erhaltenen wirk- 
lichen Wert tibereinstimmt. 
Es sei hier beilaufig bemerkt, dass die aus der Auflésungstheorie fiir lineare 


eee eee er 


Gleichungssysteme und Integralgleichungen zweiter Art hervorgehenden Neu- — 
mannschen Reihen fast immer schlecht konvergieren oder divergieren. Nun — 
verhalt sich aber eine Neumannsche Reihe fast immer ungefahr wie eine geo- — 


metrische, so dass man sie sehr wohl mit Hilfe des QD-Algorithmus sum- 
mieren kann (und zwar auch im Falle der Divergenz). Diesbeziigliche Versuche 


bei partiellen Differentialgleichungen, die allerdings erst im kleinen MaBstab ~ 


durchgefiihrt wurden, verliefen durchaus positiv. 


4. Auflésung von algebraischen Gleichungen 


Es seien die Nullstellen eines Polynoms N(z) = 2" + --- zu bestimmen. Es 
erscheint zweckmassig, zunachst eine rationale Funktion /(z) = N,(z)/N(z) zu 


konstruieren, deren Zahler irgendein Polynom vom Grade 7 — 1 ist, dann diese 
co 


in eine Potenzreihe ) (s,./z*+1) zu entwickeln und schliesslich durch Anwendung 
0 


des QD-Algorithmus auf die s, die Pole von /(z) zu bestimmen. Der QD-Algo- 


rithmus liefert diese nach I, Satz 3, als lim g“. 
a—>Co 


Leider erweist sich nun dieses Vorgehen in vielen Fallen als undurchfiihrbar, 
weil die Berechnung der «héheren» Quotienten q°, g®”, g¢?, ... fiir grosse yv in 
der Regel mit einem erheblichen Genauigkeitsverlust verbunden ist, wie das 
folgende Beispiel zeigt: 
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Beispiel: Bestimmung der Pole von 


Za 3 eee 1 ao ibe 7 
23 — 9 32_ 8 4 2 e babsagrer zB at 


‘Das zugehorige QD-Schema lautet in kompromierter Form, wenn mit sechs 
wesentlichen Stellen gerechnet wird: 


hee alam aE 9 By pens sag 
i 
at 
12). | = 225 
9,75 —0,231487 
117 0,05342 — 1,238 26 
9,803 42 — 1,52318 0,716621 
1147 — 0,008 30 0,582573 0,158641 
| 9,79512 — 0,932 307 0,292689 
11235 | 0,000 79 — 0,182893 
| 9,79591 —1,11599 
110057 | — 0,00009 
| 9,79582 
1078099 © 


Die Werte q{ konvergieren in der Tat gegen eine Nullstelle des Polynoms 
2° — 922 8z-+ 2, aber infolge der kleinen und damit relativ ungenauen 
Werte e® werden alle weiter rechts stehenden Gréssen unzuverlassig. Insbe- 
sondere wird auch ef) nicht 0, wie es auf Grund der Theorie (I, Satz 2) zu 
erwarten ware. 


5. Die progressive Form des QD-Algorithmus 


Die wenig ermutigenden Resultate des § 4 lassen sich wesentlich verbessern, 
wenn man die Formeln (4) in I dazu beniitzt, um aus einer Schragzeile des 
QD-Schemas (komprimierte Form) die nachst tiefer liegende zu berechnen. Zu 
diesem Zweck sind nur die Formeln (I, 4) umzustellen und dann unter Beach- 
tung von et) = 0 nacheinander fiir o = 1, 2, 3, ... anzuwenden: 


g* A) eee g”) we el”) sal et) f 


oO 
2) g) (4) 


(vr 1) __ o +z¢6+1 


an en . 
Wenn also einmal die in der obersten Schragzeile stehenden Elemente q{) und 
e) bekannt sind, kann man vermége der Formel (4) das ganze QD-Schema von 
oben nach unten aufbauen (statt von links nach rechts). Die progressive Form 
des QD-Algorithmus ist natiirlich besonders dann vorteilhaft, wenn die s, die 
Koeffizienten einer rationalen Funktion sind. Nach I, Satz 2, ist dann das 
QD-Schema nach rechts beschrankt (das heisst e®) = 0 fiir ein gewisses 7 und 
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alle v), so dass die Berechnung einer neuen Schragzeile nach (4) nur endlich 
viele Rechenoperationen erfordert. 

Schematische Darstellung dieses Vorgangs (die Pfeile geben an, welche 
Gréssen zur Berechnung eines neuen Elements herangezogen werden miissen) : 


gf? /o) 
7 g® 
@ 
O) | yids eo 
=) (0) 
e 
| 8 1 er | oc en-4 (0 
(1) te) 
a A | Pe | : | 
eo Seay ! ; 
ie yee i. 
ef) | 1326 In | 
= ef?) = 1 
Tn 


6. Auflésung algebraischer Gleichungen mit Hilfe des progressiven © 
QD-Algorithmus 


Wir wenden uns wieder der Aufgabe zu, die algebraische Gleichung N(z)=0 — 
durch Bestimmung der Pole von f(z) = N,(z)/N(z) aufzulésen. Der progressive 
QD-Algorithmus ist das geeignete Werkzeug fiir diesen Zweck, weil er in natiir- 
licher Weise und ohne wesentlichen Genauigkeitsverlust die g) fiir grosses 
zu berechnen erlaubt?). Auf diese Weise kann man die Nullstellen von N(z) 
beliebig genau approximieren, soweit sie nicht gleiche Absolutbetrage besitzen. 
Ein Problem bietet lediglich die Beschaffung der Anfangswerte g!°) und e). 
Man findet diese nach (2) offenbar durch Entwicklung der rationalen Funktion 
N,(z)/N(z) in einen S-Kettenbruch mit Hilfe des Euklidischen Algorithmus: 

Ausgehend von den Polynomen 


Ny (2) = 214+ --- und NF(z2) = Nz) =2"4+--- 


konstruieren wir fiir k = 1, 2,3, ... nacheinander die Polynome N,(z) und N#(z) 

1) Nachtrag bei der Kovvektur. Von Herrn Dr. A. HOUSEHOLDER, Oak Ridge, Tenn., werde 
ich darauf aufmerksam gemacht, dass A1rken fiir seine unter [1] behandelte Erweiterung der 
Bernoullischen Methode ebenfalls eine progressive Form angibt. Vgl. A. C. A1rKeNn, Further 
Numerical Studies in Algebraic Equations and Matrices, Proc. Royal Soc. Edinburgh 5/, 80-90 


(1931). Uber die sich daraus fiir den QD-Algorithmus ergebenden Folgerungen wird spater be- 
richtet werden. 
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yom Grade n — k: 

| 2 N,(2) — N&i(2) 
qT 

Ne (2) — Nyle) 


ex, 


ic 


N*(2) = 


(eae re | ce 
| 5 


J ale) = (k= 1,2,...%—1). 
Dabei werden die q, und e; eindeutig bestimmt, wenn man verlangt, dass die 
Koeffizienten der héchsten Potenzen in allen diesen Polynomen N, und N# 1 
‘sein sollen. 

Da man nachtraglich leicht feststellt, dass der endliche Kettenbruch 


on 

die Funktion /(z) darstellt, sind diese q,, e, gerade die gesuchten Anfangswerte 
gi) e), wie ein Vergleich mit (2) zeigt. 
Beispiel: Bestimmung der Pole von 


Be ae Sige 1 
EE aoe 


Wir erhalten nach (5), wenn man mit sechs wesentlichen Stellen rechnet (aus- 
gehend von N,* = 2° — 9 22 —- 8z+ 2; N,= 224+ 32+ 1): 


a Nj = 1222+ 92-2; also q, = 12 N¥ = 22 + 0,75 z — 0,166667 
Ns = — 2,25 s —1,16667 also e, = — 2,25 N, = 2+ 0,518519 

q,2N$ = — 0,231481 z+ 0,166667 also g, = —0,231481 N¥ = z—0,720003 

@,.N, = — 1,23852 é, = —1,23852 N,=1 

G3. N$ = 0,720003 q3 = 0,720003 N*¥ =1. 


Damit kann das QD-Schema nach (4) aufgebaut werden. (Die s, werden fiir 
diesen Zweck gar nicht gebraucht und sind darum weggelassen) : 


tt 4 ey” a” ey qs” (5) 
Y 12 
— 2,25 
9,75 —0,231481 
0,0534187 —1,23852 
9,803 42 —1,52342 0,720003 
— 0,008 3011 0,585 353 0 
9,79512 — 0,929766 0,134650 
0,000 7880 — 0,0847716 0 
9,79591 — 1,015 33 0,219422 
— 0,000 0816 0,018 3199 0 
9,795 83 — 0,996 928 0,201 102 
: + 0,000 008 3 — 0,003 695 5 0 
‘ — 1,00063 0,204 798 
. 0,000 7563 0 


0,204012 
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Das Bou zeigt, dass die numerische Instabilitat beseitigt ist; die drei Ko- 
lonnen 9°”, q%, 7%) konvergieren offensichtlich gegen die gesuchten Pole 


J, = 5 +23 w 9,79583; Ag——1 und ”25=5—V23 ~ 0,204168. 


Fiir die Falle, wo die Folgen der g®) nicht geniigend gut konvergieren, sind 
konvergenzverbessernde Massnahmen vorgesehen.. Sie beruhen darauf, dass 
man auf die unendlichen Reihen in der aus (4) folgenden Beziehung 


fore) co | 
A, = tim gf) SqFt DT a7 = DF er (6) 
Tee!) “n=v x=v41 d 


die Methoden von § 3 anwendet. 
Wenn die zu lésende Gleichung konjugiert komplexe Wurzeln hat, so zeigt | 
sich dies durch ein Oszillieren gewisser g-Werte an, genauer: Wenn 


A, | = | Ag| = ers > |An—1| == |A;,| > |Anst| = reat 


so konvergieren zwar die Kolonnen g® fiir o + k —1,k nach wie vor gegen 
Wurzeln der Gleichung, aber die Quotienten g), und gq?) oszillieren, und ins- — 
besondere konvergiert e®), fiir »-—>oo nicht gegen Null. Trotzdem finden © 
wir in I, Satz 4, und or eeconaere in Folgerung 2 ein Mittel zur Bestimmung 1 
der konjugiert komplexen Wurzeln A,_, und A; (man setze dort m = k — 2 und © 
Hie R)\ 


Man bilde die Polynome 
BY (2) = 2? — [gD + Pp ]z t+ 90) as, (6a) 
die fiir vy > oo gegen (z — Ay_y) (2 — A,) Ronvergieren. 
Beispiel: Auflésung der Gleichung 24 — 4 23 + 5 22 — 42+ 1 durch Bestim- ~ 


mung der Pole von 29/(z4 — 4 23 + 5 z2 4241). Der S-Kettenbruch dieser — 
Funktion wurde nach (5) bestimmt und lautet: j 


1] Ad gnli2sl 045). orn ey ater 

Zz aa as 2 | , t= ida ager tii, Ses ’ oa 

ASS evare Arr haar ie bl x od rae 
0,0292207|  0,392.857 | 


ae te 1 


| 
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Damit lasst sich das QD-Schema aufbauen (in komprimierter Form und ohne 
die s®)): 


5 q® e) q’? el?) q ew) q 
4 
—1,25 
PAST KS) 0,45 
— 0,204545 — 0,977778 
2,545 45 — 0,323 232 1,41414 
0,025 9740 4,27778 — 0,029 2207 
2,57143 3,92857 — 2,892 86 0,392 857 
0,0396825 — 3,15000 0,003 968 2 
2,61111 0,738 889 0,261111 0,388 889 
0,011 2293 —1,11316 0,0059102 
2,622 34 — 0,385496 1,38018 0,382979 
— 0,0016508 3,985 39 0,001 6400 
| 3,601 54 — 2,603 57 0,381 339 
— 2,881 06 — 0,000 2402 
Ay 0,277 730 0,381579 
0 — 0,000 3300 
Se ee y 0,381 909 
komplexes Wurzelpaar ? 0) 
Ay 


Da g und q? nicht zu konvergieren scheinen, bilden wir nach (6a) die Poly- 
nome #” (mit k = 3): 

£0) z* — 1,09091 z + 0,636364 , 

pO) — 22 — 1,03571 z + 0,935065 , 

be) yee + 1,02579 , 

p) — z2 — 0,994.68 z + 1,01980, 

pit) <a\2® — 0,99797 z + 1,00367 , 


I 


Diese konvergieren fiir » > oo offenbar gegen z? — z + 1, welches in der Tat 
ein Teiler von z4—422+522—42+41 ist; die zugehdrigen Wurzeln sind 


fog = (1 + i V3)/2. 
7. Quadratische Konvergenz des QD-Algorithmus 


Man kann die in § 5 behandelte progressive Form des QD-Algorithmus auch 
als eine Methode zur Berechnung der S-Kettenbriiche fiir die in I, § 5, defi- 


nierten Funktionen 
iej=)>) Sete feet oe ets.) 
0 


sus dem S-Kettenbruch von f,(z) = /(z) auffassen. Wir wollen nun auch noch 
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die J-Kettenbriiche dieser Funktionen in den Kreis unserer Betrachtungen 
ziehen, wobei wir der Einfachheit halber den S- (bzw. J-) Kettenbruch der 
Funktion /,(z) mit S, (bzw. J,) bezeichnen. Auf diese Weise wird die Konver- 
genz des QD-Algorithmus wesentlich beschleunigt werden kénnen. 

Man kann namlich die Formeln (4), die dem progressiven QD-Algorithmus 
zugrunde liegen, offenbar in zwei Formelgruppen aufspalten: 


ale) = 9b + of | 


fiir o = 1, 2;...n, mit &) = 0; (7) 
Berd = go, | 
6 Lol, 6 
gt) = at) — een | . 
way OD firo—1,2,...,.”, mité@*?=P'1=0. (8)§ 
Coe ger) } 
o 


; 
Wie die Ausfiihrungen von I, §5, zeigen’), sind die hier auftretenden af *), BY** — 
tatsadchlich die Elemente des /-Kettenbruchs /,,,, der demnach als Zwischen- © 
glied der Berechnung von S,,, aus S, erscheint. : 

Man kann nun mit Vorteil bentitzen, dass bei einer als /-Kettenbruch gege- — 

benen Funktion die Verschiebung des Koordinatennullpunktes eine triviale 
Operation ist, man muss nur zu jedem «, dieselbe Konstante addieren, denn es ~ 
ist ja 

B cdeet sire hn Be he ake pee eal eee 


2 Oy le — a, | z =a, 


wenn z=2+4, he = %, +t. 

Es ist nun zweckmissig, im Zwischenschritt, das heisst unmittelbar nach 
der Berechnung der «*! und B°*”, den Nullpunkt der z-Ebene in die Nahe | 
eines Pols von f(z) zu verschieben (wir bezeichnen diesen Pol fortan mit /,). 
Unter geeigneten Umstanden wird der Absolutbetrag dieses Pols im neuen Ko- 
ordinatensystem wesentlich kleiner sein als bei den andern Polen. Das wird 
nach den Ausfithrungen des Beweises von Satz 3 in der Arbeit I zur Folge 
haben, dass die e®), im weiteren Verlauf der Rechnung wesentlich rascher 
gegen 0 konvergieren. Das hat wiederum zur Folge, dass man bessere Nahe- 
rungswerte fiir A, erhalt, die man fiir weitere Koordinatenverschiebungen be- 
nutzen kann, so dass die Konvergenz immer schneller wird. 

Fihrt man bei jedem Schritt der obigen Rechenvorschrift eine solche Koor- 
dinatenverschiebung aus und wird jeweils g®) als Naherungswert fiir A, im 
betreffenden Koordinatensystem gewahlt, so ergibt sich (nachdem man die 


1) Man beachte insbesondere den Kettenbruch 


iC ee [! 4 ci Se aes See geal =) 2 oy fl) 


s—qi-e{ = |z—g3—-e3 | z—-qh a 2 Zz 
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®) und p°*” wieder eliminiert hat) die folgende Rechenvorschrift fiir die 

Berechnung der Nullstellen eines Polynoms N(z): 

1. Zuerst wird wie in § 6 eine rationale Funktion N,(z)/N(z) konstruiert und 
_ sogleich in einen S-Kettenbruch entwickelt. 

2. Ausgehend von den Elementen q®), e dieses Kettenbruchs und mit f, = 0, 

berechne man nacheinander fiir y = 0, 1, 2,... 


1=44+¢, 
Z (y+1 Ng 
| co eae ln te Gaal Baa 0) 
(@) 40) were 
Bey <2 4o+1 fo ree ‘ 
a eal) mit ¢ "= ¢, =0. 


Dann konvergieren die e®) gegen 0 und die Gréssen ¢, + g”) gegen die gesuchten 
Nullstellen von N(z), and 2War konvergieren e®), und ¢, + ¢g® quadratisch 
gegen 0 bzw. A,. 
_ Sobald e®), hinreichend klein geworden ist, lasst man gq’) und e“), weg und 
fiihrt die Pe hear orschrilt mit ~ — 1 an Stelle von weiter. Auf diese Weise 
erhalt man auch /,,_, mit quadratischer Konvergenz usw. Die Praxis zeigt, dass 
die Berechnung von 4, relativ lange dauert, wahrend dann die tibrigen A, 
rasch nacheinander erhalten werden. Wenn auch allfallige komplexe Wurzeln 
mit quadratischer Konvergenz erhalten werden sollen, so muss man auch kom- 
plexe Koordinatenverschiebungen ausfihren. 

Numerisches Beispiel. Wir bestimmen wie in § 6 die Nullstellen von 

3 __ 9 z2 — 8z + 2 als Pole der rationalen Funktion 


ze+3¢+1 
fh Sie ae Y 


Der S-Kettenbruch lautet hier nach § 6: 


1 12> | 225 0, 231481 — iz oe 0,720003 | 
oe se a ee eC 
Damit ergibt die Rechenvorschrift (9) bei sechsstelliger Rechnung: 
) a ey ro ey qs ly 
) 12 2 — 0,231481 — 1,23852 0,720003 0 
L 9,03000 0,057678 — 2,24768 0,396 735 — 0,396735 0,720003 
2 9,48441 — 0,013 669 — 1,44054 0,109 264 — 0,109 264 0,323 268 
3 9,58005 0,002055 — 1,22407 0,009 753 — 0,009753 0,214004 
I 9,591 81 — 0,000 262 — 1,204 30 0,000079 — 0,000079 0,204 251 
; 9,591 63 0,000033 —1,20417 0 0 0,204172 
10,7958 — 0,000 004 0,000 004 — 0,999 998 
10,7958 0) 0) — 0,999995 


ro ae ee Ar 


0 9,79580_ 
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Die unterstrichenen Werte sind die gesuchten Nullstellen. Die Abweichungen 
von den wahren Werten sind auf Aufrundungsfehler zuriickzufiihren, insbe 
sondere auch auf solche wihrend der Berechnung des S-Kettenbruches (mar 
beachte zum Beispiel den ungenauen Anfangswert 93 = 0,720003 an Stelle 
von 0,720000). 


Das Literaturverzeichnis befindet sich am Schlusse dcr Arbeit I (H. Rutis- 
HAUSER, Der Quotienten-Differenzen-Algorithmus, ZAMP 5, 233 [1954]). 


Summary 


The present paper gives several applications of the quotient-difference 
(= QD)-algorithm [see ZAMP 5, 233 (1954)]. These include the summation of 
infinite series with the aid of continued fractions and the solution of algebraic 
equations. For the latter, a slight modification of the QD-algorithm is developed 
in § 5. 


(Eingegangen: 18. September 1953.) 
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Ein Impulsgenerator mit Sekundaremissionsroéhren 


Von Ernst BaLDINGER und Marc NICcoLet, Basel?) 


Der Entwurf dieses Apparates wurde durch eine Arbeit von F. H. WELLS?) ~ 
angeregt. Der Apparat erzeugt Rechteckimpulse einstellbarer Dauer sowie Span- — 
nungen zur Zeitablenkung und zur Hellsteuerung eines Oszillographen und wird 
zur Kontrolle und zur Messung der charakteristischen Eigenschaften von Breit- 
bandverstarkern, Laufzeitketten usw. verwendet. Die grundsatzliche Anordnung 
ist dem Blockdiagramm (Figur 1) zu entnehmen. 

Ein Kippgenerator G liefert eine Sdgezahnspannung Z und einen Rechteck- 
impuls P, der — tiber die Kathodenstufe A — zur Hellsteuerung der Kathoden- 
strahlréhre dient. Der Impuls P steuert ebenfalls die Schaltung J, in welcher der 
Testimpuls S erzeugt wird. Das RC-Glied V gestattet, diesen Testimpuls inner- 
halb der Zeit 7, zu verschieben. Die Ausgangsstufe A liefert sowohl positive als 
auch negative Impulse einstellbarer Grésse. 

Die Impulse P und S werden mit einer Kippschaltung erzeugt, die als Grund- 
element eine Sekundaremissionspentode verwendet. Ihr Vorteil besteht darin, 
dass infolge der gleichzeitigen Anwesenheit von positiven und negativen Signalen 
alle impulsmassig gesteuerten Rdhren mit Ruhestrom Null betrieben werden 
konnen. Infolgedessen lassen sich iiberall grosse Momentanwerte der Stréme und 
somit kurze Anstiegzeiten der Impulse erreichen. Das Prinzip der verwendeten 
Kippschaltung ist in Figur 2 dargestellt. Im Gegensatz zu iiblichen Réhren fiihrt 


1) Physikalisches Institut der Universitat Basel, Abteilung fiir angewandte Physik. 
) F.H.Wetts, Nucleonics 10, Nr. 4, 28 (1952). I. A.D. Lewis und F. H. We tts, Milli- 
microsecond Pulse Techniques (Pergamon Press, London 1954), 
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ie€ Anode mehr Strom als die Kathode, weshalb die Schaltung zu Kippvorgangen 
ahig ist. Die Réhre kann, einmal in den leitenden Zustand ubergegangen, auf 
erschiedene Arten wieder gesperrt werden: 

. Durch die Entladung des Kondensators C,.; 

- durch ein negatives Gittersignal, wie dies beim Generator G der Fall ist; 

. durch eine in die Anode eingefiihrte, am Ende kurzgeschlossene Laufzeitkette, 
deren Reflexsignal die Kathode wieder tiber das Gitterpotential hebt und so 
die Impulsdauer bestimmt. Nach diesem Prinzip arbeitet der Generator IE 
der den Testimpuls S erzeugt. 

Um einer Uberlastung der Sekundaremissionsréhren vorzubeugen, ist das Ver- 
haltnis 7, : 7, in Figur 1 mit 1 : 100 fest vorgeschrieben. Man erreicht dies durch 


Figur 1 Figur 2 
Blockdiagramm. Prinzip einer Kippschaltung mit Sekundaremissionsréhre. 


jie in Figur 3 wiedergegebene Anordnung. Fiihrt namlich die Réhre 7 Strom, so 
wird der Kondensator C’ durch die Pentode 2 angenahert linear entladen. Dieser 
Vorgang hért auf, wenn die Riickfiihrdiode 5 zu leiten beginnt, wodurch beide 
RGhren gesperrt werden. Uber die Triode 3 wird nun C’ mit einem rund hundert- 
mal kleineren Strom wieder aufgeladen. Sobald der Ansprechpegel der Schmitt- 
Schaltung 41) iiberschritten wird, fiihrt diese dem Gitter der EF P60 (Roéhre 7) 
sin positives Signal zu, wodurch die Entladung des Kondensators C’ von neuem 
seginnt. Die fiir den ganzen Zyklus benotigte Zeit hangt vom Werte von C’ ab. 
Jm eine méglichst lineare Entladung des Kondensators C’ zu gewahrleisten, ist 
ss zweckmAassig, dem Bremsgitter der Rohre 2 ein gegentiber der Kathode posi- 
ives Potential von ~ 50 V zu erteilen. Diese Massnahme verbessert die Ja-Ua- 
Kennlinie der Pentode und erhGht insbesondere bei kleinen Anodenspannungen 
hren Innenwiderstand?). 

Das Signal P wird nun einerseits der Kathodenstufe /y (Hellsteuerung), an- 
lererseits einem RC-Glied zugefiihrt und so aus dem Rechteckimpuls eine expo- 
ventiell ansteigende Kurve erzeugt. Ihr Schnittpunkt mit dem Ansprechpegel 
ron J bestimmt den Zeitpunkt des Einsatzes des Testimpulses in /. Der Dreh- 
condensator C gestattet, diesen Punkt innerhalb von 7, zu verschieben. 

Durch Anderung der Gliederzahl der Laufzeitkette in der Kippschaltung / 
isst sich die Dauer des Testimpulses zwischen zirka 10~* und 127-1077 s ver- 


1) Vgl. zum Beispiel W. C. E-more und H. Sanps, Electronics (McGraw-Hill, New York, 1949), 


. 100. ‘ 
2) Vgl. zum Beispiel K. R. SpancenBERG, Vacuum Tubes (Mc Graw-Hill, New York, 1948), 


}. 272, 
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62 EERO OB 

2 : 
3 12AX7 S zur : 
Z 12AT7 {lL Endstufe A — 
5 tf 6H6 
7 6AL5 | 


R,=1kQ, R, = 8002, R,= 4902, Ry=2kQ, R=5—160kQ, Zy = 330 Q, 
Cy =0,5uF, Cy=100pF—oo, Cy=200pF, C= 5—550pF, 
C’ = Schaltkapazitat — 6000 pF 
Figur 3 


Schaltschema der Teile G, / und J von Figur 1. 


ee a ee 


a 
10* s 
Figur 4 


Testimpuls von 1077 s Dauer. 


andern. Figur 4 zeigt den kiirzesten Impuls von 10-7 s Dauer, aufgenommen auf 
dem Philips-Oszillographen Typ GM5653/02 in Kombination mit der raschesten 
Zeitablenkung des Impulsgenerators. 

Testimpuls und Zeitablenkspannung werden direkt den Ablenkplatten der 
Kathodenstrahlréhre zugefiihrt. Die Ausgangsstufe A in Figur 1 — ebenfalls eine 
EF P60 — liefert im giinstigsten Falle einen Impulsstrom von = 0,28 A oder eine 
maximale Impulsspannung von & 400 V. Zum Philips-Oszillographen sei noch 
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bemerkt, dass die Verdrahtung zu den Ablenkplatten der Braunschen Rohre ab- 
geandert werden muss, da sonst stérende, hochfrequente Schwingungen auftreten?) 
__ Zur Beobachtung rasch ablaufender Vorgange wird der Abfall (Z,, Figur 1), 
ftir langsam ablaufende Vorgange der Anstieg (72, Figur 1) der Sagezahnspan- 
mung verwendet. Die Repetitionsfrequenz ist einstellbar zwischen 5-104 und 


10% Hz, was einer Variation der Schreibgeschwindigkeit von 5: 107 bis 104 cm/s 
entspricht. 


Summary 
j ; 4S 4 5 , | 
A pulse generator employing secondary emission valves is described. The 


generator produces positive and negative output pulses with adjustable pulse 
duration from 10~7 to 10-*s and a rise time of about 1,3: 10-8 s. In addition to 
these test-pulses, signals for the beam brightening and the sweep of a cathode-ray 


oscillograph are generated, permitting the study of high-speed transients with 
an ordinary oscilloscope. 


(Eingegangen: 13. Juli 1954.) 


On Rotational Viscous Flow Through a Tube’) 
By PETER HENRIcI, Washington, D.C.) 


1. Introduction 


In a recent paper [3]*), L. Cottatz and H. GOrTLER consider the stationary, 
slightly rotational, axisymmetric flow of a viscous fluid through a straight cir- 
cular tube at moderate or high Reynolds numbers. It is shown that after some 
simplification this problem can be reduced to the following eigenvalue problem 
for a function G(s) proportional to the angular velocity component: 


SG te Ss Ge G =wAs* (1 s4\)Gs, US 5s Sil, (1a) 
Ci GA e: (1b) 


With the aid of the eigenfunctions of this problem the rotation in the tube can 
be computed for any given initial rotation. Furthermore, the n-th eigenvalue of 
(1) measures the exponential decay of the corresponding component of rotation 
in the direction of the flow. 

The first five eigenvalues and eigenfunctions of (1) are computed in [3] by 
numerical methods, and a crude approximation for the large eigenvalues is given. 
In the present note we supplement these results by asymptotic expressions for 
the large eigenvalues and the corresponding eigenfunctions. It turns out that in 

1) Die tiefste Eigenresonanz der abgeanderten Zuleitungen (inklusive der Kapazitat der Ablenk- 
platten) liegt bei 80 MHz. Der aperiodische Grenzwiderstand betragt ~ 150 Q. Mit einer Verdrah- 
tungskapazitat von ~ 20 pF entspricht dies einer Anstiegszeit von ~ 0,5 + 1078s. 

2) This paper was prepared under a National Bureau of Standards contract with American 
University and was supported by the Office of Naval Research. 

3) American University. 

4) Numbers in brackets refer to References, page 516. 
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particular for the large eigenvalues we are thus able to get much better approxi 
mations than those obtained in [3]. We also indicate an alternative method fo 
the computation of the low eigenvalues.— The physical conclusions of [3] are no 
affected by our results. e 


2. Analytical Solution of the Eigenvalue Problem 


In the following we put 
A= ALO BP (2 


It is easy to verify that the differential equation (1a) and the first of the bound 
ary conditions (1b) are satisfied by the function 


G(s) eBieie 7 es 1F,(1 — B; 2; 4B s*), 


I 


gre M5 y2(4 pst 


where the symbols ,/, and M denote the confluent hypergeometric series and 
the Whittaker function!) respectively. The eigenvalues of the problem can then 
be computed by (2) from the nonvanishing zeros of the function 


f(B) = Mg,1/2(4 B) . 


Taking for granted from [3] that the eigenvalues are all negative, we can confine 
our search for these zeros to real and, since?) 


H(—8) = —#(8), . 


to positive values of . 


3. Integral Representation for G(s) 


The last factor on the right of (3) is, as is easily seen, equal to the coefficient — 
of z~1 in the Laurent expansion at z = oo of the function : 


1 ei bss (1 “ia =)’. : 
B ; 


which in the z-plane cut between = 0 and z= 1 is univalued and regular 
except at » = oo. By the residue theorem we thus have ; 


e72bs 


: Das: (* 2, 1\B E 
c= AEE fete hae al 


where the integral is to be extended along an arbitrary simple contour encircling 
the mentioned cut clockwise. Setting . 


g(e) = 258 4542 — log(1— —) (7)) 
we can also write : 
oe 
Diels eee 1) E 
Cts) = / e PB) gy | (8) 


1) In Wurrraker and Watson’s original notation; see E. T. WuirrAKER and G.N. Warson, 
A Course of Modern Analysis, 4th ed. (Cambridge University Press, London, 1927), p. 337. 
*) Modern Analysis, p. 338. 


rf 
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4. Large Eigenvalues 


In order to obtain information about the eigenvalues, we consider first the 
case s = 1. In this case the integrand in (8) possesses the two (only seemingly 
coinciding) saddle points of the second order , 


if : 
och as +70 (9) 
on the upper and the lower edge of the cut. At these points, 
g= Fina, ois Oy g”=0, gv + 0. 


Integrating along the appropriate lines of steepest descent originating from 2,» 
by integrating from 0 to co with respect to the new variable ; 


4: ie ae 


2 ge Oa. Pe 
« “ ae 


j= 


and defining the coefficients b,, (m= 1, 2,...) by 


wen foe) 
a em rh Din (3 i lau 1) | 
m=1 (10) 
3 3 
(b= 1, B= Sn be ae, ea = 0) | 


we obtain by a well-known lemma due to Watson®) the following asymptotic 
expansion for the function /(f)%), valid for | arg 6| S «—0 with arbitrary 6>0: 


i\p) ~ — = O's cosa (6+ = os a r(3) ey 
a ee | re 
es ts) Ge)" coe (° os *) t aH ("+ i 
\3 
Bet ne. 
; cos (P— 3) ae 70 r(2} Eames els 


Using an argument similar to one put forward by Watson‘), we can derive from 
(11) the first terms of an asymptotic expansion for the large zeros of f(B)*). The 


1) The branch of €(t) which is real for positive ¢ is to be considered here. 
2) G.N. Warson, 4 Treatise on the Theorie of Bessel Functions, 2nd ed. (Cambridge University 


Press, London, 1944), p. 236. 
3) The first terms of this expansion can also be found in [1], p. 109. They are there in a form 


which is not very well suited to our purpose. 


4) Bessel Functions, p. 503. : 
5) A similar expansion for the large zeros of Mg (4 B) has been given by LAuwERIER [4], also 


in connection with a problem of viscous flow, 


ZAMP V/33 
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ama. Js ta 
result is as follows: 
i! 1 i 1 \ - 10/3 F ’ 
Ba Wate + Ag ag Ao [(™ es 3 | (12) 
(n+ =) (n+ %) 


where x is a larger integer and 


1/3 7 


= Gh we 


We do not prove that £, approximates just the n-th zero of f(f); this should be 
clear, however, from the numerical results given below. Formula (12) should be 
contrasted with the approximation given in [3] for the large eigenvalues, which 
in the present notation amounts to 


bam nN. 


A= 


5. Qualitative Behaviour of the Eigenfunctions 


In order to discuss the asymptotic behaviour of the eigenfunctions, we assume 
in(8) (0) s <= 7 Puttine 
S=— cos? (0 < Os 5), 


we find that the saddle points lie now at 


1 bits 


— 


41, fee 


bo 


and are simple. Without further amplification we state the result of the saddle- 
point integration, which is in accordance with [1], p. 105: 


| 
: 
| 
. 
: 


pe ke he ee Wy? whe sce 1 
G(s) ~ rE sin ( (x — 28+ sin2 0) — +). | ~ 


(e<osF-e, B oo}. | 


This formula suggests that the relative maxima of |G(s)| are smallest near 


oe = u , i.e. s = 21/7 and increase from there on towards the endpoints of the 
considered interval. In this direction we prove the following 


Theorem. Let B > 1/8 and 


aa 3 i oe | 
Sj ee ie 4 16 pe ’ ie SON 


Then the sequence of the relative maxima of |G(s)| decreases in (s_, So) and in- 
creases im (So, S;) . : 


Proof. Equation (1a) may be written in the selfadjoint form 


GaGa = 
where Pao sgt we 


Ce ie 


~ »* 
ie. w% 
4 


Since g(s) > 0 in (sii, S,) and since s g(s) is increasing in (s_,, So) and decreasing 
in (So, $;), the assertion follows from a theorem due to SonINE and Porya), — It 
will be noted that 
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nh 


1 
$4 ~ 7 wand eae 


4B 
6. Low Eigenvalues 


_ The eigenvalues of low order are computed in [3] as zeros of a series which is 
defined by a three-term recurrence relation. It follows from (3) that the correspond- 
ing values of f may also be determined as the zeros of the confluent hypergeo- 
metric series 


1k = 6; 25 4 By =D y's (14) 
n=0 


the terms of which may be computed by the two-term relation 


n— Bp 
Pe ms a ig a ep n(n ra) Yn-1- 


7. Numerical Results 


The following is a table of the first few values of f,, as computed from the 
values of 4, given in [3], computed by the asymptotic expansion (12) and com- 
puted exactly from the zeros of (14). The calculation of (14), from which the zeros 
were determined by the regula falsi, was performed on the National Bureau of 
Standards Eastern Automatic Computer (SEAC) with the aid of the standard 
subroutine (working with double precision) for the hypergeometric function 
described in [2]*). 


By 
n 
[3] (12) (14) 

1 1.156025 1-155760 1-156025 

2 2-1619 2-161833 2-161876 

3 3-163 3+ 163742 Se WOStoD 
| 4 4-146 4-164635 4-164640 
| 5 (308 4 5-165140 5-165143 


We have not computed an estimate for the remainder in (12) and are therefore 
unable to give bounds for the error involved in using the asymptotic formula*). The 
slose agreement already between the first exact eigenvalues and the values com- 
puted by the asymptotic formula should however be noted. 


1) See SzEG6 [5], p. 161, footnote. 
2) It is a pleasure to acknowledge the able assistance of Mr. W. F. Canixt in the performance 


of these computations. ; 
3) In this respect also the argument leading to the estimate for the error of J, on p. 106 of [3] 


s not quite complete, because there is no indication how the lower bound —/, of —A, was derived. 
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Zusammenfassung 


Fiir die Eigenwerte und’Eigenfunktionen eines Eigenwertproblems, auf das © 
L. Cottatz und H. G6rRTLER im Zusammenhang mit einer kiirzlichen Unter- — 
suchung [3] tiber Rohrstr6mung mit schwachem Drall gefiihrt worden sind, wer- 
den asymptotische Naherungsformeln angegeben, die insbesondere die héhern 
Eigenwerte mit grésserer Genauigkeit als in [3] zu berechnen gestatten. Auch zur 
Ermittlung der niedrigen Eigenwerte wird ein einfaches Verfahren vorgeschlagen. 
Die Ergebnisse werden numerisch mit denen von [3] verglichen. 


(Received: June 10, 1954.) 
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Internationaler Mathematikerkongress in Amsterdam 
September 1954 


Der in Abstaénden von vier Jahren tagende Internationale Kongress fiir 
Mathematik war in diesem Jahr von den hollandischen Mathematikern — speziell 
vom Prasidenten, Prof. SCHOUTEN, und dem Sekretar, Prof. KoxsMa — in hervor- 
ragender Weise organisiert worden. Er verreinigte die fiihrenden Mathematiker 
der meisten europaischen Lander, der nord- und siidamerikanischen Staaten, der 
USSR., Australiens und der Staaten des ferneren Osten. In sieben parallel ge- 
fiihrten Sektionen (Algebra und Arithmetik, Analysis, Geometrie, Statistik und 
Wahrscheinlichkeitsrechnung, Angewandte und numerische Mathematik, Grund- 
lagen und Logik, Unterricht) wurden die wissenschaftlichen Mitteilungen entge- 
gengenommen. Ausserdem waren verschiedene Gelehrte (worunter zwei Schwei- 
zer; ECKMANN und STIEFEL) von der Kongressleitung zu langeren Berichten iiber 
ihre Spezialgebiete eingeladen worden. An der Eréffnungssitzung sprach Prof. 
VON NEUMANN, Princeton, tiber Ungeliste Probleme dey Mathematik, und Prof. 
WEYL iibergab die sogenannte «Field-Medaille» an die beiden Mathematiker 
Koperra und SERRE. Diese Medaille wird an junge Mathematiker verliehen, die 
in den verflossenen vier Jahren besonders zur Entwicklung der Mathematik bei- 
getragen haben. 

In der Sektion fiir angewandte und numerische Mathematik traten vor allem 
die stiirmischen Umwertungen zutage, welche der Einsatz der neueren mathema- - 
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tischen Instrumente (elektronische Rechenmaschinen usw.) mit sich gebracht 
hat. Die einstiindigen Hauptvortrage dieser Sektion waren: 


Prof. VAN DANT7z1G, Holland, iiber die mathematischen Probleme, welche durch die 
Uberschwemmungskatastrophen des Jahres 1953 in Holland angeregt wurden. 

GOLDSTEIN: On Some Methods of Approximation in Fluid Mechanics. 

JESSEN: Some Aspect of the Theory of Almost Periodic Functions. 

NEYMAN: Current Problems of Mathematical Statistics. 

STIEFEL: Recent Developments in Relaxation Techniques. 

TitcHMARSH: Eigenfunction Problems Arising from Differential Equations. 

‘Yosipa: Semi-group Theory and the Integration Problem of Diffusion Equations. 


Halbstiindige Vortrage hielten: 
Coiiatz: Fehlermassprinzipien in dey praktischen Analysis. 
FicHERA: Methods of Funktional Linear Analysis in Mathematical Physics. 
HEsTENES: Hilbert Space Methods in Calculus of Variations and Numerical Analysis. 
‘KAMPE DE FEriEtT: Problémes mathématiques posés par la mécanique statistique 
de la turbulence. 
Rewiicn: Halbbeschrinkte Differentialoperatoren hoherer Ordnung. 
STOCKER: Recent Progress in the Theory of Surface Waves. 
Es ist eine grosse Ehre fiir die Schweiz, dass Prof. Dr. H. Hopr fiir die nichste 
Amtsperiode als Prasident der Internationalen Mathematischen Union gewahlt 
wurde. E. STIEFEL 
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Einftihrung in die Physik, Band 1: Mechanik, Hydvomechanitk und Thermo- 
dynamik. Von P. FRAUENFELDER und P. HuBER (Ernst-Reinhardt-Verlag AG., 
Basel 1951). 492 S., 391 Abb.; sFr. 18.50. 

Das vorliegende Buch ist nach dem Vorwort als Erganzung zu Vorlesungen 
iiber Experimentalphysik gedacht. Dies erlaubt, grésstenteils auf die Beschrei- 
bung von Apparaten und Versuchen zu verzichten, und gibt dafiir Gelegenheit 
zur Konzentration und zu einem systematischen Stoffaufbau. Ein solcher bedingt 
stets eine eindimensionale Ordnung der weitraumigen Stoffgebiete der Physik, 
wofiir sich eine in jeder Hinsicht ideale Lésung wohl iiberhaupt nicht finden lasst. 

Im Buche von FRAUENFELDER und HUBER werden in einem einleitenden 
Kapitel zunachst die Grundeinheiten von Lange, Zeit und Masse mit den Dimen- 
sionsbegriffen eingefiihrt. Daran schliesst sich eine kurze Darstellung der Vektor- 
algebra und der Fehlertheorie. Die Mechanik beginnt darauf mit der Statik der 
starren Kérper unter Voranstellung der Newtonschen Prinzipien. Es folgen die 
Dynamik des Massenpunktes und der Systeme, die Elastizitatslehre und die 
Hydromechanik. Der zweite Teil des Buches ist der Thermodynamik gewidmet. 

Beim Lesen spiirt man die grosse Erfahrung der Verfasser in der Vermittlung 
von Lehrstoff. In wohltuendem Gegensatz zu manchen andern Lehrbiichern wird 
auf Punkte, die erfahrungsgemass den Studierenden Schwierigkeiten bieten, mit 
besonderer Sorgfalt eingegangen. Gut geraten in dieser Hinsicht sind vor allem 
einzelne Kapitel der Thermodynamik und der kinetischen Theorie der Warme. 

Wenn man auch gegen Einzelheiten Vorbehalte anbringen kénnte — so etwa 
gegen die von den Verfassern getroffene Wahl der mechanischen Grundbegriffe —, 
so bedeutet das keine Einschrankung der Anerkennung des hohen Wertes des 


Buches und der Klarheit und Exaktheit seiner Darstellung. Je mehr man s 
~ darin vertieft, um so mehr iiberrascht die Fiille des Stoffes, die wohl selten 
‘dieser konzentrierten und doch immer klar verstandlichen Form geboten wird. 
Besonders wertvoll diirften fiir den mit dem Stoff weniger Vertrauten die zahl- 


_. reichen ergiinzenden Zahlenbeispiele sein, bei denen je nach dem Stoffgebiete alle 


drei heute iiblichen Ma®Bsysteme (CGS-, Giorgi- und technisches System) zur 
Anwendung gelangen. a 
Die Sprache des Buches ist die knappe, sachliche Sprache des Physikers. Die 
- sehr zahlreichen, trefflich gezeichneten, meist schematischen Abbildungen stellen 
eine wertvolle Erganzung dar. Drucktechnisch wirkt der Text dank der Verwen- 
dung verschiedener Typen sehr iibersichtlich. Zu bedauern ist, dass fiir die Gros— 
 senbezeichnung nicht konsequent kursive Buchstaben Verwendung fanden. Auch» 
_ beanspruchen leider die zur Charakterisierung der Vektoren beniitzten Pfeile 
etwas viel Platz und bringen wegen des dadurch bedingten wechselnden Zeilen- 
_ abstandes eine gewisse Unruhe in das Textbild. 
Gemass dem Vorwort ist das Buch fiir naturwissenschaftlich interessierte Stu- 
dierende gedacht. Es vermag aber auch bedeutend héhere Anspriiche zu befrie- 
_digen und diirfte jedem, der sich mit physikalischen Fragen beschaftigt, ein wert- 
volles Hilfsmittel fiir ein vertieftes Eindringen sein. Man wird daher das Erschei- 
nen des zweiten Bandes mit besonderer Spannung erwarten. W. Hardmeier_ 
Advances in Electronics. Vol. IV, Edited by L. MArton (Academic Press, © 
New York, 1952). X + 344 pp., 118 figs.; $7.80. “gy 
Uber Ziel und Anlage der von L. MARTON herausgegebenen Sammlung «Ad- 
vances in Electronics» haben wir bereits bei der Besprechung der friiheren Bande 
dieser Biicherreihe berichtet. Der vorliegende vierte Band umfasst die folgenden — 
Beitrage (welche diesmal bis auf einen aus amerikanischen Forschungsstatten | 
stammen): H. S. W. Massey, Electron Scattering in Solids; G. A. Morton, The ~ 
Scintillation Counter; A. VAN DER ZIEL, Fluctuation Phenomena; C. V. L. SmitTH, ~ 
Electronic Digital Computers; J. S. Donat, Jr., Modulation of Continwous-Wave 
Magnetrons; W. E. Fromm, The Magnetic Airborne Detector; M. G. PAWLEY and ~ 
W.E. Triest, Multichannel Radio Telemetering. S| 
Was von den friiheren Banden Lobendes gesagt werden konnte, kann auch — 
fir Band IV wiederholt werden. Jedenfalls darf auch dieser Band auf das Inter- — 
esse eines weiten Leserkreises zahlen, bildet doch jeder einzelne der Beitrage eine © 
ausgezeichnete Ubersicht iiber das von ihm behandelte Gebiet. A.A. Rusterholz — 


Sampling Techniques. By W.C. CocHran (J. Wiley & Sons, New York, — 
1953), 330 pp.; $6.50. 

Dieses Buch gibt eine gute Darstellung der modernen Stichprobentheorie — 
und -technik, wie sie hauptsachlich in den letzten zwanzig Jahren, unter wesent- 
licher Mitwirkung des Autors dieses Buches, entwickelt wurde. Es behandelt vor 
allem die technischen und praktischen Fragen, die sich bei der Anwendung stati- _ 
stischer Stichproben stellen, wie zum Beispiel: Wie gross soll die Stichprobe 
gewahlt werden, und welche Schliisse kénnen aus ihr gezogen werden? Welche 
Sicherheit besitzen diese Schliisse ? Die Theorie wird stets an Hand praktischer — 
Beispiele demonstriert; ferner enthalt das Buch eine Fiille von Aufgaben. 

Der Leser muss die wichtigsten Verteilungen, ihre zugehérigen Begriffe und — 
statistischen Teste kennen, wenn er das Buch verstehen will. Es gibt allen denen 
eine vorztigliche Anleitung, die auf méglichst rationelle Art Statistiken anlegen | 
und deren Resultat verwerten wollen. W. Saxer ~ 
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Neuerscheinung 


Schaltungstheorie und Messtechnik | 
des Dezimeter- und Zentimeter- 
Wellengebietes 


von Dr. Albert Weissfloch, Technischer Direktor der Sté. Techn. d’Application et de 
Recherche Electronique, Montrouge (Seine). 


(Sammlung Lehrbiicher und Monographien aus dem Gebiet der exakten Wissen- 
schaften, Physikalische Reihe, Band 1.) 


Inhalt: Allgemeine kreisgeometrische Vierpoltheorie — Schaltungen mit homogenen 
Leitungen (Doppel- und Wellenleiter) und Messtechnik — Zusammengesetzte Schal- 
tungsteile, (2)-Pole mit m > 2 (Verzweigungen, Richtungskoppler) und Frequenz- 
abhangigkeit der Transformationseigenschaften (Resonatoren, Filter) - Anpassung. 


Der erste Abschnitt gibt eine ausfiihrliche Darstellung der Geometrie der linear 
gebrochenen Funktionen, die, angewandt auf die Vierpoltheorie, ein sehr niitzliches 
mathematisches Werkzeug zur Lésung von Hochstfrequenzfragen bietet. Die wei- 
teren Abschnitte behandeln dann speziellere Schaltungsteile, mit denen sich der’ 
Ingenieur befasst. 
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